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Abstract: We study the charge transport properties of fields confined to a (2+l)-dimensional 
, defect coupled to (3+l)-dimensional super- Yang-Mills at large-iV c and strong coupling, using 

AdS/CFT techniques applied to linear response theory. The dual system is described by N t 
qq | probe D5- or D7-branes in the gravitational background of iV c black D3-branes. Surprisingly, 

the transport properties of both defect CFT's are essentially identical - even though the 
D7-brane construction breaks all super symmetries. We find that the system possesses a 
conduction threshold given by the wave-number of the perturbation and that the charge 
transport arises from a quasiparticle spectrum which is consistent with an intuitive picture 
where the defect acquires a finite width. We also examine finite-A modifications arising from 
higher derivative interactions in the probe brane action. 

X 

Keywords: AdS/CFT correspondence, Gauge/gravity correspondence. 



Contents 



1. Introduction 2 

2. Defect branes 3 

2.1 D5-branes 4 

2.2 D7 probes 8 

3. Correlators 9 

3.1 Spectral functions 10 

3.1.1 Transverse correlator 13 

3.1.2 Longitudinal correlator 14 

3.2 Temperature and / dependence 16 

3.2.1 Effective temperature 17 

3.2.2 Resonances on the defect 19 

4. Quasinormal modes and quasiparticles 22 

4.1 Finding the Ansatz 23 

4.2 Quasiparticles from the collisionless regime 25 

4.3 Poles in the hydrodynamic regime 29 

5. Electromagnetic duality and perturbative corrections 31 

6. Hall conductivity 36 

7. Discussion 38 

A. Diffusion constant on the defect 41 

B. Corrections to D5-brane action 42 

C. Hyperbolic tangent potential 46 

C.l Finding the spectral curves 46 

C.2 Exact pole structure 48 



- 1 - 



1. Introduction 

The AdS/CFT correspondence provides a powerful framework for the study of strongly cou- 
pled gauge theories [1, 2, 3]. While the gauge theories that are currently amenable to such 
holographic analysis are typically very different from real world QCD, this duality may still 
give us some insight into the strongly coupled quark-gluon plasma (sQGP) produced in recent 
experiments at RHIC [4]. Matching AdS/CFT results with experimental data here antici- 
pates that the sQGP can be described by an effective CFT which is in the same universality 
class as the strongly coupled gauge theories studied holographically. Similar considerations 
have recently motivated exploring the possible application of the AdS/CFT correspondence 
to the study of (2 + l)-dimensional condensed matter systems [5]. A variety of holographic 
models displaying interesting properties, including superfluidity, superconductivity and Hall 
conductivity, have now been studied [6, 7, 8]. Further interesting models of various types 
of nonrelativistic CFT's have also been constructed [9]. One advantageous aspect of the 
AdS/CFT correspondence is the "uniformity" of the calculations, i.e., a single set of calcu- 
lations describes the system in different disparate regimes (to — ► versus T — > 0). This can 
be contrasted with more conventional field theory analysis of conformal systems in (2+1)- 
dimensions [10]. One might also approach these explorations of possible "AdS/Condensed 
Matter" correspondences with some optimism because of the rich variety of experimental 
systems which might exhibit conformal behaviour and also the vast array of four-dimensional 
AdS vacua in the string theory landscape [11]. 

In the present paper, we use the AdS/CFT correspondence to study the charge trans- 
port properties of fields confined to a (2+l)-dimensional defect coupled to (3+l)-dimensional 
super- Yang-Mills. The defect CFT is realized by inserting N f probe D5- or D7-branes into 
the background of a black D3-brane. In either construction, we are also able to introduce 
an additional internal flux, in which case the defect in the CFT separates regions where the 
rank of the SYM gauge group is different. The defect CFT constructed with the D5-branes 
is certainly well known [12, 13, 14]. Certain aspects of the D7-brane construction have also 
been studied previously [15, 16] but we should note that the internal flux introduced here is 
essential to remove an instability that would otherwise appear in this construction. We will 
use linear response theory to study the conductivity on the defect at finite frequency and 
temperature and at finite wave-number, i.e., the conductivity of an anisotropic current. 

An overview of the paper is as follows: In section 2, we review the holographic frame- 
work, in particular the embedding of the probe D-branes in the AdS5 x S 5 background. In 
section 3, we obtain the basic results for the spectral functions, starting with a review of the 
methodology and then the computation of the transverse and longitudinal conductivities in 
section 3.1. Here we also comment on the agreement with the diffusion-dominated conduc- 
tivity in the hydrodynamic regime, i.e., in the regime at small frequency and wave- number 
uj,q «T. This is followed by a discussion of the collisionless, (q T), regime using analyti- 
cal approximations in section 3.2, for both the insulating case (at small fequencies) and the 
optical regime (at large frequencies). Using those results, we study the spectrum of quasi- 
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normal modes in section 4. In section 5, we examine the effect of stringy corrections to the 
gauge theory on the probe branes, which describes the behaviour of the dual currents on the 
defect. In particular, electromagnetic duality is lost when these a'-corrections are included, 
which has interesting implications for the conductivity as strong but finite 't Hooft coupling. 
Finally, we consider the computation of a topological Hall conductivity in section 6. Section 
7 closes with some discussion and observations about our results. Some details of our analysis 
are relegated to appendices: In appendix A, we calculate the diffusion constant for charge 
transport. Appendix B presents some details of the analysis including certain ol corrections 
in the D5-brane worldvolume action. In appendix C, we do an analytical study of a slightly 
simplified model of the defect, which gives further qualitative insight, and aids the numerical 
computation of the quasinormal modes. 

2. Defect branes 

The AdS/CFT correspondence is most studied and best understood as the duality between 
type lib string theory on AdS§ x S 5 and M = 4 super- Yang- Mills theory with U(N C ) gauge 
group. In this context, all fields in the SYM theory transform in the adjoint representation of 
the gauge group. One approach to introducing matter fields transforming in the fundamental 
representation is to insert probe D7-branes into the supergravity background [17]. However, 
this approach can also be used to construct a defect field theory, where the fundamental fields 
are only supported on a subspace within the four-dimensional spacetime of the gauge theory. 
In particular, we will consider constructing a (2 + l)-dimensional defect by inserting N f Dp- 
branes, with three dimensions parallel to the SYM directions and p — 3 directions wrapped 
on the S 5 . In the following, we work with both probe D5- and D7-branes. If we consider the 
supergravity background as the throat geometry of iV c D3-branes, our defect constructions 
are described by the following array: 



The D5-brane construction is supersymmetric and the dual field theory is now the SYM 
gauge theory coupled to N f fundamental hypermultiplets, which are confined to a (2+1)- 
dimensional defect. Note that the supersymmetry has been reduced from J\f = 4 to N = 2 by 
the introduction of the defect. In the D7-brane case, we have lost supersymmetry altogether 
and the defect supports N f flavours of fermions, again in the fundamental representation [15]. 
One should worry that the lack of supersymmetry in the latter case will manifest itself with 
the appearance of instabilities. However, we will explicitly show below in section 2.2 that this 
problem can be avoided. In the limit N f <C N c , the D5- and D7-branes may be treated as 
probes in the supergravity background, i.e., we may ignore the gravitational back-reaction of 
the branes. 



background 
probe 



0123456789 

D3 x x x x 

D5 x x x xxx 

D7 xxx xxxxx 



(2.1) 
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As we commented above, a similar holographic framework has been used extensively 
to study the properties of the N = 2 gauge theory constructed with parallel D7- and D3- 
branes, i.e., the fundamental fields propagate in the full four-dimensional spacetime - e.g., 
see [18, 19, 20]. If a mass M q is introduced for the hypermultiplets, it was found that the 
scale M iun ~ M q /\/A plays a special role in this theory. First, the "mesons" , bound states of a 
fundamental and an anti-fundamental field, are deeply bound with their spectrum of masses 
characterized by M fun [21]. Next at a temperature T ~ M fun , the system undergoes a phase 
transition characterized by the dissociation of the mesonic bound states [20]. The analogous 
results can be verified for the defect theories considered here. That is, the meson spectrum is 
characterized by the same mass scale M iun [22, 23] and these states are completely dissociated 
in a phase transition at T ~ M fun [24]. However, these results are tangential to the present 
study, as we will only consider the conformal regime with M q = 0. 

Common to both of our constructions is the supergravity background dual to N = 4 
SYM at finite temperature. This background is a planar black hole in AdS*,, corresponding 
to the decoupling limit of iV c black D3-branes [25]: 

ds 2 = ^ (-h(r)dt 2 + dx 2 + dy 2 + dz 2 ) + ^ (J?L + r 2 dn 2 ^ , c[% z = - £ (2.2) 

where h(r) = 1 — r 4 /r 4 . The gauge theory directions correspond the coordinates {t, x,y, z}. 
The radius of curvature L is defined in terms of the string coupling constant g s and the string 
length scale £ B as L 4 = 4irg B N c £f. The holographic dictionary relates the Yang- Mills and 
string coupling constants as g 2 M = 4ng s and so we may write L 4 = X£f where A = g 2 M N c 
is the 't Hooft coupling. As usual, we work in the supergravity approximation, ignoring 
the effects of string loops or higher derivative terms suppressed by powers of £ s (except in 
section 5 and appendix B). Hence, we are working in the limit where both N c , A — > oo. The 
background (2.2) contains an event horizon at r = ro- The temperature of the SYM theory 
is then equivalent to the Hawking temperature: 

T-£. (2.3) 

2.1 D5-branes 

Introducing D5-branes as in (2.1) was the original application of probe branes for the holo- 
graphic construction of a defect CFT - e.g., see [12, 13]. The worldvolume action which will 
determine the embedding of the probe D5-branes has the usual Dirac-Born-Infeld (DBI) and 
Wess-Zumino (WZ) terms: 

I 5 = -N { T 5 J dPa^-det (P[G\ + 2ir£ 2 F) + N f T 5 f A 2< 2 F . (2.4) 

Implicitly we have assumed that the N f D5-branes are all coincident. Hence, in principle, their 
worldvolume supports a U (N f ) gauge theory, however, implicitly above and in the following, 
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we only consider the gauge field in the diagonal U(l) of this U(N { ). We choose coordinates 
on the five-sphere in (2.2) such that 

dfi§ = <# 2 + cos 2 $ (d9 2 + sin 2 6 d(f) + sin 2 $ dfL 2 . . (2.5) 

The D5-branes wrap the two-sphere parameterized by {#,</>} above, fill three of the gauge 
theory directions {t, x, y} and extend in the radial direction r. We also introduce a flux of 
the worldvolume gauge field on the two-sphere: 

ify = 2^sin0. (2.6) 

One may verify that this flux corresponds to dissolving q D3-branes into the worldvolume 
of the N { D5-branes along the {t, x, y, r} directions, since the branes with flux sources 

C(4) 

through the WZ term in (2.4). 

Now in general, the D5-brane embedding would be specified by giving its profile in both 
the angular direction ip{r) and the D3-brane direction z(r). These embeddings all have 
translational symmetry in the {i, x, y}-space, as well as invariance under SO (3) rotations on 
the internal two-sphere. In the following, we consider only the embeddings with tp = 0, i.e., 
where the D5-brane wraps a maximal two-sphere in the internal space. One can easily verify 
this choice corresponds to a solution of the worldvolume embedding equations. This choice 
also corresponds to setting the mass of the fundamental fields to zero, i.e., M q = 0, and so 
as we will describe below, this choice also ensures that the dual field theory with the defect 
remains conformal. 

Hence in our analysis, we must determine the profile z(r). The induced metric on the 
D5-branes is now described by 

r 2 / I 2 r 2 \ 

ds 2 = (-h(r)dt 2 + dx 2 + dy 2 ) + ( + J2 d ^\ dr2 + ^ ( d02 + sin2 9 ^ ■ ( 2 - 7 ) 

We can integrate over the two-sphere directions to produce a factor of 

d 2 ayJdet S 2 (P[G] + 2n^F) = 4ir (iV 2 L 4 + ir 2 tfq 2 ) 1/2 = 4vriV f L 2 ^1 + / 2 (2.8) 



i 



is 2 

where 



f=<±-^L± (2 9) 

in the DBI part of the action (2.4). The full D5-brane action then becomes 

I 5 = -4TrN f T 5 ^/T+ f 2 J d 3 xdrr 2 (l + ^h(r)d r z 2 y - 4irN f T 5 -j^ J d?xdrr A d r z. 

(2.10) 

To simplify the analysis, we introduce the following coordinates: 

u=—, X = T2 Z - ( 2 -ll) 

r L A 
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With this new notation, h(u) = 1 — u 4 and so the horizon is now at u = 1 while the asymptotic 
region is reached when u — > 0. The worldvolume action can now be written as: 



du 



h = -47rrgiV f T 5 / d 3 z ^ + P (l + %)x' 2 ) 5 + / x' 



(2.12) 



where %' = <9 u x- This expression is independent of %, such that the variation with respect to 
X' yields a constant of motion: 



1 



u 1 



h(u)x' 



(1 + h(u) X 



mV 2 



+ / 



c. 



(2.13) 



To avoid singular behaviour at the horizon, we need to fix the integration constant to be 
C = f. In this case, (2.13) yields 



X 



f 



2„4 



(2.14) 



Given this expression, the profile x( u ) can De expressed in terms of an incomplete elliptic 
integral. However, in the following, it will sufficient to have a closed form expression for 
We illustrate some typical profiles in figure 1. 

In terms of our original coordinates, 
we have 



o . o 



dz 
dr 



L 2 f 



2 12 



(2.15) 




Figure 1: The brane profile \ f° r various values 
of / at finite (solid) and zero (dashed) temperature. 
(Note the T — profiles correspond to x = — f u -) 



Here we may consider the supersymmet- 
ric limit with r$ = 0, in which case (2.15) 
simplifies to z = —L 2 f/r. In this case, one 
can confirm that the induced metric (2.7) 
corresponds to AdS4 x S 2 with an AdS ra- 
dius of curvature of LAdS 4 = L \fl + f 2 
[12]. Hence the system inherits SO(2,3) 
symmetry from the AdS4 geometry, which 
reflects the fact that the dual field theory 

remains conformal in the presence of the defect. This conformal invar iance can also be shown 
directly by an analysis of the field theory [26]. Subsequently, the construction of the fully 
back-reacted geometries corresponding to the D5-branes embedded in AdS$ x S 5 demonstrate 
that the preservation of the SO(2,3) symmetry is a fully nonperturbative result [14]. 

One may note that with the supersymmetric profile, z = —L 2 f/r, there are an additional 
q = y/XNff/ir D3-branes stretching from z = to — oo, assuming / > 0. Hence if one were 
to include back-reaction, the asymptotic five-form flux would be shifted from iV c — > N c + q 
units on this side of the space. The same will apply at a finite temperature. Even though 
the brane falls through the horizon at a finite distance in this case, continuity at the horizon 
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dictates that the background will carry N c + q units of flux out to z = — oo. In either case, 
the natural interpretation is that the dual CFT has a U(N C + q) gauge group in the region 
z < 0, while the gauge group remains U(N C ) for z > 0. 

It is interesting to pursue the interpretation of the above brane configuration in the dual 
CFT further. A detailed AdS/CFT dictionary has been developed for this defect system 
[13, 26]. In particular, one finds that the defect lagrangian contains potential source terms 
for the adjoint scalars in the SYM theory [13, 27]. The D5-brane carrying flux / corresponds 
to producing a noncommutative configuration of adjoint scalars in a U(q) subgroup of the 
U(N c + q) in the z < region [22]. In fact, in this super symmetric configuration, the profile of 
the D5-branes can be precisely matched to the scalar profile using noncommutative geometry 
[28]: r 2 = (2< 2 ) 2 ^Tr(<I> 2 ) where ^Tr(d> 2 ) = ^ 

As the D5-brane wraps a maximal two-sphere inside the 5 5 , one might worry about the 
stability of this configuration. Indeed the worldvolume field, corresponding to fluctuations in 
the angle ip, is found to be a tachyon with [12] 

^--ppW-ZE^-nL- (2 ' 16) 

However, the last inequality indicates that the ■0-mode satisfies the Breitenlohner-Freedman 
bound [29] in the asymptotically AdS^ geometry induced on the D5-brane worldvolume. 
Hence this field does not in fact produce an instability. 

Another concern may arise in considering the intersection of the D5-branes with the event 
horizon in (2.2). There we note that 

_f_ 



*1«=i = — >=S ( 2 - 17 ) 



or in terms of original coordinates 



dz 
dr 



L 2 f , s 

(2.18) 



=r 



Since the D5-brane enters the event horizon at an angle, one might worry that the induced 
geometry is singular [30]. However, one can verify that this intuition is incorrect and that 
in fact, the D5-brane geometry remains smooth as it crosses the horizon. Hence the induced 
metric (2.7) describes a smooth 'black hole' geometry on the D5-brane worldvolume. A related 
question is: what is the surface gravity or the temperature of the induced horizon? It is a 
simple exercise to show that the relevant temperature matches that of the bulk geometry, 
i.e., that given in (2.3). Of course, this reflects the fact that the defect and bulk fields will 
be in thermal equilibrium, as expected. 

We address one other potential concern related to the internal flux (2.6). Throughout 
the paper, we will be considering finite values of /, typically of O(l). Hence according to 
(2.9), we are introducing q ~ 0{yfX) D3-branes and so one might worry about whether it is 
reasonable to consider the probe brane limit, i.e., to ignore the gravitational back-reaction of 
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the branes. Of course, this is not a problem since the overall tension of the D5-branes is not 
significantly modified by the flux, as can be seen from (2.10). The essential point is that the 
D3-branes are distributed on the D5-branes over the internal two-sphere which has an area 
of order L 2 ~ y/\ and so the density of D3-branes remains small, i.e., the density is 0(f). 

2.2 D7 probes 

The case of D7 probe branes is similar to the previous section with D5-branes. The main 
difference lies in the internal part of the geometry. In particular, the D7-branes wrap a(n 
equatorial) four-sphere in the internal S 5 . As before, we consider D3-branes dissolved into 
the probe branes. In the present case, the D7-branes source the three-brane charge through 
in the appropriate term in the WZ action: ^(27r£ 2 ) 2 T 7 J A F A F. Hence, considering a 
stack of N f coincident D7-branes with a U (N { ) gauge symmetry, we introduce a nonvanishing 
second Chern class on the internal four-sphere: q 7 = ^ § s4 , TrF A F. 

The D7-branes are fixed to wrap a maximal four-sphere while the embedding in the AdSs 
is described by z = z(r). The induced metric on the D7-branes becomes 

ds 2 = ^ (-h(r)dt 2 + dx 2 + dy 2 ) + (^^y + J^' 2 ) dr 2 + L 2 dflj . (2.19) 

Since the present configuration contains a nontrivial nonabelian gauge field, the worldvolume 
action requires a nonabelian extension of the DBI action [31] 

I 7 = -T 7 J d 8 a STrV -det (P[G] + 2n£ 2 F) + ^T 7 (2irl 2 ) 2 J Tr (c^ AFAfj . (2.20) 

This action uses the proposal of a maximally symmetric gauge trace, denoted by 'STY' [32]. 
To be precise, the trace includes a symmetric average over all orderings of F a i, - and implicitly 
any appearances of the nonabelian scalars as well [33] but the latter will not be relevant in the 
present analysis. This prescription correctly agrees with the string action to fourth order in 
the field strength [32] but is known to miss certain commutator terms which begin to appear 
at sixth order [34]. However, the contribution of such terms is typically suppressed by factors 
of 1/N ( and so they can be safely neglected for sufficiently large N t [28, 35]. 

As before, we integrate over the internal space in the DBI action. Here the internal S* 4 
carries an nonabelian gauge field giving the instanton number q 7 . This configuration was 
extensively studied in [35] and hence using their results, we find 



Js 



d A n STr y/det S 4 (P[G] + 2ir£ 2 F) 

s 4 

= jfyn Vfli^ 8 + \L\2v£ 2 ) 2 F ah F^ + ^(2^2)4(^^06^)2 

= ^-( Nt L' + ^ 2 £t\q 7 \) . (2.21) 

In the latter, we use (anti-)self-duality for the instanton configuration: F a i = (— )\e a bcdF cd 
for q 7 > (q 7 < 0). Implicitly, we are also assuming that the instanton number is uniform on 
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the four-sphere, which limits q 7 < N f (N 2 — l)/6 [35]. Substituting (2.21) and the embedding 
(2.19) into (2.20), the action for the background configuration becomes 



h = -^L\l + \Q\)T 7 J Sadr^l + ^Kr) - ^L 4 QT 7 J Sodr^z* , (2.22) 

where we defined for convenience Q = 67r 2 ^-^ = ^f-|^- Now, the computations analogous 
to those in section 2.1 yield an identical embedding x{ u ) as in (2-14) but the constant / is 
replaced by 

h = ; Q • (2-23) 

The microscopic interpretation of the D7-brane configuration in the dual CFT is not as clear 
in the present case. However, as before, the gauge group in the region z < will be enhanced 
to U(N C + q 7 ) assuming q 7 > 0. There should be source terms on the defect which excite a 
noncommutative configuration of the adjoint scalars in the transverse space. The latter can 
be interpreted in terms of noncommutative geometry as giving the profile of the D7-branes, 
at least to leading order in l/iV f [35]. 

An important difference between the present case and that in the previous section with 
D5-branes is that in the mass of the tachyonic mode ip corresponding to the S 4 part of the 
D7-branes "slipping off" the maximal 5 4 in the internal space. A simple calculation reveals 
that 

4 _ 4 \Q\ + l (2 24) 

^ ~ L*(\Q\ + 1) " Li dS4 2|Q| + 1 " [2 - M) 

Recall that the BF bound requires m 2 > — jft — [29] and hence is only satisfied for \Q\ > |. 

Hence one can trust the results in the following sections for the D7-branes only for /| > 49/32 
and we might think of the internal flux on the S 4 as creating some pressure that stabilizes 
the size of the S 4 . However, we should caution the reader that what we have shown is that 
the most obvious instability is removed for sufficiently large Q. While suggestive, this does 
not prove the D7-brane configuration is absolutely stable. 

Beyond this crucial difference, the analysis of these two systems (i.e., defects constructed 
with D5- or D7-branes) is completely the same. Hence in the following, we focus on the 
first case of D5-branes and only comment on differences in coefficients that may arise for 
D7-branes where appropriate. 



3. Correlators 



In this section, we obtain examine various correlators of the currents dual to the worldvolume 
gauge field A^. First we review the basic form of the correlators below, following [5]. Then 
we numerically compute the spectral functions in 3.1 and then examine the dependence of 
the correlators on the temperature and the flux / in 3.2. 



- 9 - 



In the following, we use holographic techniques to calculate the retarded Green's function 
for a conserved current J^(x) on the defect. The defect degrees of freedom form a (2+1)- 
dimensional CFT which restricts the form of the correlators: 

C^{x -y) = -i 9(x° - y°) ( [J^x), My)} > , (3.1) 

where translation invariance is assumed. The correlator can be Fourier transformed to C^ v {p) 
with p^ = (to, k)} Now current conservation and rotational invariance (full Lorentz invariance 
is lost with T > 0) restrict the form of the Fourier transform of this correlator to be [5] 

C^ip) = Pj, Il T (p) + P^ v U L (p) . (3.2) 

where the transverse and longitudinal projectors can be written as 

pT _ X. . _ kik 3 pT _ n 

If we take into account that the conformal dimension of the current J^(x) is 2, the components 
U T ' L in (3.2) take the form: 

II T ' L (p) = y^K T ' L (lo/T, k/T). (3.4) 

In the limit of T = 0, we have H T (p) = Ii L {p) = II(j>) and recover the Lorentz invariant 
correlator 

c^(p)= (^-^?) n(p). (3.5) 

In order to produce physical observables, and to interpret our results from a condensed matter 
point of view, we will calculate the conductivity from the Kubo formula 

&ij = . (3.6) 

3.1 Spectral functions 

In this section, we compute spectral functions for excitations of fundamental fields on the 
defect by studying fluctuations of the world volume fields on the D5-brane probes. In particu- 
lar, we focus on correlators of the the world volume vector A^, which is dual to the conserved 
current J M corresponding to the diagonal U(l) of the global flavour symmetry on the defect. 
The worldvolume gauge field gives rise to several types of modes, one of which is a vector with 
respect to the Lorentz group in the (2+l)-dimensional defect. These modes are characterized 
as having only -Ao,i,2 nonzero while the components on the internal two-sphere are vanishing 
[23]. Further the radial component A r can consistently be set to zero because we only study 
modes which are constant on the internal space [23]. 

x We work with the mostly positive signature so that rj^ — diag( — 1, +1, +1). 
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While the full action for the gauge fields on the D5-branes receives contributions from 
both the Dirac-Born-Infeld (DBI) action plus a Wess-Zumino term, since our gauge field 
fluctuations have vanishing radial and S 2 components, only the DBI portion of the action is 
relevant in determining their dynamics. Since we only study linearized fluctuations about the 
background, the gauge field action is only needed to quadratic order, which is simply 

(2< 2 ) 2 r2 _ 



L gauge 



-4ttL 2 y/1 + f 2 N { T 5 J d 3 a dr^g~ ■ 



1 



/ d A odr^gF l . (3.7) 



4ft 2 



Here we have integrated over the internal S 2 as in (2.8) and use g^ to denote the induced 
metric (2.7) in the AdSs directions. Above, we also defined the effective gauge coupling for 
the four-dimensional Maxwell field: 



9\ 

For the D7 case, this becomes 

1 _32vr 4 iV f/)4r4 



9l 



t s L\l + \Q\)T 7 = (l + \Q\) 



N t N c 
3vr 2 



(3.8) 



(3.9) 



Note that the gauge field action (3.7) corresponds to the standard Maxwell action in a four- 
dimensional curved spacetime. Hence, these gauge fluctuations will exhibit electromagnetic 
duality, which was shown to play and interesting role in the physics of the conformal field 
theory in [5]. We will explore this point further in section 5. Of course, Maxwell's equations 
follow as 

d a (y=g~F ab )=ti. (3.10) 

Using these equations of motion, the Maxwell action (3.7) becomes a total derivative 
and following the standard prescription, we obtain the desired correlator from the resulting 
boundary term. To proceed, let us first give the explicit metric on the brane, 



T 2 

ds 2 = K 



r l ( - h{u)dt 2 + dx 2 + dy 2 ) + H (1 + h(u) X ' 2 ) 



(3.11) 



in terms of the dimensionless radial coordinate u, as given in (2.11). Then the action (3.7) 
becomes 



L gauge 



'2g 2 J 



d 3 a du d a 



~gA b F ab ] = ~ J d 3 a [V^g aa g uu A a d u A a ] u u Zl 



1 f f 1 

= ~2gjh\ dZ ° (l + h(u) X l2 ) l l 2 ^ u)AxduAx + h (u)A y d u A y - A t 8 u A t 



U-t-l 

(3.12) 



The usual AdS/CFT prescription tells us that we will only need the contribution at the 
asymptotic boundary u — > [36]. Following [37], we take the Fourier transform of the gauge 
field, 

—e ik -°A„(k,u), (3.13) 
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to write the boundary action as 
1 1 f 



L gauge 



2gjl (2vr)3 „ 
2g 2 V (2vr)3 J 



lU-t-l 



d 3 k [^-gg^g uu A^u,-k)d u A,(u,k)\ u _ V! 

(1 + h(u) X ' 2 r 1/2 (Hu)A x (u, -k)d u A x (u, k) 



d 3 k 



+h(u)A y (u, -k)d u A y (u, k) - A t (u, -k)d u A t (u, fc))]J£j 

with a single sum of \x being implicit in the first line. 

Looking at the asymptotic behaviour of the fields, we write 



(3.14) 



A^(k,u) = A^{k) 



A^(k,u) 
A^(k,u ) ' 



(3.15) 



where no is a UV regulator and it is understood that eventually the limit no — > will be 
taken. We can then derive the flux factor for, say, A y by taking variations with respect to 
A y0 [36]: 



T — 

•'yy ~ 2 



h(u) 



Ay(u,-k)d u A y (u,k) 



where e = ^2 



9a 



a + h{u)x' 2 y/ 2 Ay(u ,-k)A y (u ,k) (,U(,) 

we will show later how £0 relates to the charge permittivity. The 
flux (3.16) should be conserved, i.e., be independent of the radius n. The usual AdS/CFT pre- 
scription tells us to evaluate it at the asymptotic boundary, while applying infalling boundary 
conditions at the horizon (n = 1), to find the retarded Green's function (3.1) for the current 
J M in the defect CFT [36]: 



a 



yy 



-2K 



yy 



£0 



h{u) 



A y (u,-k)d u A y (u,k) 



_ (1 + /i(u)x' 2 ) 1/2 A/K —k)Ay(uo, k) 



ll,«()^0 



£0 



d u A y (u, k) 



(3.17) 



The other correlators C^ v follow in general by rewriting (3.15) as A^(k, n) = A l/ o(k)M u /1 (u, k) 
[5] and making the variation SA p SA - . In our case the t, t and x,x correlators are given by 
(3.17) with the indices appropriately replaced: 



Cx 



d u A x (u,k) 



and C, 



it 



d u A t (u, k) 



(3.18) 



In order to evaluate the spectral function, we must solve the equations of motion (3.10). 
It is convenient to introduce dimensionless coordinates by rescaling the defect coordinates as 



t = ^t 
L 2 



x 



L 2 X > 



ro 

y = T2 y. 



(3.19) 



Without loss of generality, we also assume the fluctuations only carry momentum in the x 
direction, i.e., k^ = (Co, q, 0) — note that, e.g., Co = L 2 jr§ to = lo/ttT. We note that given the 
Fourier transform (3.13), the vector potentials vary as e lkl±x>i in the gauge theory directions. 
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Now the explicit equations of motion simplify to 



b = 


u : 









(3.20) 


b = 


- 1 : 





* 2H 1 h 


(q 2 A~ t + tiqAx) , 


(3.21) 


b = 


x : 





~ Ax+ \h 2H) 


A 'x+tf {u 2 Az + uqAi) , 


(3.22) 


b = 


y ■ 





- A y + [h 2H 




(3.23) 



where 'prime' denotes d u and 

H(u) = 1 + h{u)x 2 ■ (3.24) 
Before proceeding further, we make the following convenient definition for the conductivities 

i 

&ij = —Cij = irTcTij . (3.25) 

UJ 

Comparing to (3.6), here we are simply dividing by the dimensionless frequency uj, rather 
than uj. 

3.1.1 Transverse correlator 

Let us look carefully at the y equation (3.23). We see firstly, that in the limit u — > 0, it 
reduces to 

= A^ + (l + f 2 )(u 2 -q 2 )Ay. (3.26) 

The solution of interest is then Ay = Ago e _l v / ( 1 +/ 2 )( w2 ~9 2 ) M ) where the sign in the exponential 
is chosen so that the solution corresponds to an infalling wave. Given this solution, if one 
now calculates the correlator with (3.17) and applies (3.6), the resulting conductivity is 



T yy ~ £ o 



y/1 ~ q 2 /Co 2 . (3.27) 



The cut in the conductivity at Co = q may be surprising and we return to this point in section 
4.2. We will refer to this simple result as the low temperature approximation, reasoning as 
follows: The result applies for large dimensionless "frequencies", i.e., |(l+/ 2 )(cl> 2 — q 2 )\ 1 ^ 2 ~> 1 
to be precise. However, recalling that e.g., Co = lo/ttT, if we fix the dimensionful quantities 
{co,q}, then eq. (3.27) should apply in the limit of very low temperatures. 

To solve for the full spectral functions, we must proceed with numerical calculations. 
First, we impose infalling boundary conditions at the horizon — recall that the time-dependence 
of the potentials is e~ tujt . If we expand about 1 — u — > + , we find an appropriate description 
of the field to be 

Ay~{\- u A r' A (1 + 0(1 -«) + ■■■ ) (3.28) 

where 

i,3 + 5/ 2 q 2 / .Co s 

rTTp + V ~ l 2. 



0=^^ + 7^— A^-^)- (3-29) 
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In order to implement the infalling boundary condition and to ensure numerical stability, we 
choose the Ansatz 

Ay = (1 - U y^ 4 e-P u f(u) , (3.30) 

and solve for CF(u), which is nonsingular at the horizon, with d u CF{u) = at u = 1. As the 
second boundary condition, we fix the asymptotic normalization: A y \ u= \ = 1. 




Figure 2: The normalized transverse conductivity Re a yy (oj)/eo at q <E {1,2,3,4} for various values 
of the flux /. The low temperature approximation (3.27) is shown as the dotted line for each q. 

Figure 2 shows Re a yy (Co)/eo for various values of q and f 2 . We see that at / = the 
spectral functions are similar to those in [5]. However, as / increases, they approach the low 
temperature limit (3.27) more closely, and show some oscillatory behaviour at Co > q. We will 
discuss aspects of this behaviour in section 3.2. 

3.1.2 Longitudinal correlator 

Now, let us consider the tt equation (3.21). It is easy to see that (3.20-3.22) are not inde- 
pendent, and we cannot produce a second order equation involving Af only. However, we can 
produce one for A'~: 

- 4" - (£(f)' + §)4' + (£§(!)' -(§)')4 + f (^)4, (3.3D 

which simplifies to 

A'r \" fh! H'\ f A'r\' H o A'r 

' ! +17-^7 H= +^ 2 ~hf)^=, (3-32) 



and hence is the same as (3.23) for A y replaced by A~/s/H. Let us set A'- = c\J~H A y with 
some constant c to be determined. Now, to find d u A^(u, k) / A^(u, k) at u = 0, we employ 
(3.21) and A'~ = c\fH{A'~ + $gAy) as in [5]. It follows then from h'\ u=0 = = H'\ u=0 that 



c = ^i +/2 (^ A to + ^ A -o) • (3.33) 

A y\u=0 
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Hence, we can read off 



d u Af(u,k) 
M u > k ) 



and 



d u At(u,k) 



from 



A', 



y/1 + P^fHAy 



A'y\ U = 



Finafly as in [5], we find 

Ctt = —£q <f/C, 



uu 



and 



{t A t + CoqA x0 ) ■ 



C X x — —£q Co 2 1 Cj 



yy • 



(3.34) 



(3.35) 



Applying (3.25), these results yield interesting relations for the corresponding conductiv- 
ities. In particular, (3.35) yields 

&xx = tt/Vyy ■ (3.36) 



We can also consider a low temperature limit as above. However, this is most easily derived 
by combining (3.27) and (3.36) to find 



(3.37) 



Again, we will return to discuss the cut appearing in the conductivity at Co = q in section 4.2; 
and the conductivity can only be found in general from numerical calculations. Some typical 
results for (the real part of) shown in figure 3. We note that again that as the flux 

/ increases, our results approach the low temperature approximation (3.37), together with 
some "oscillatory" behaviour similar to that found in the transverse case. In contrast to the 
results in the previous section, the conductivity here diverges as Co — > q, as can be anticipated 
from (3.37). 



Re<T x 



6. 
5. 
4 . 
3. 
2 . 
1 . 



T -» approximation 










f = : ' 










: f- 2 :i 










f = 5 I 










; — r hi 



















0. 



2. 



3. 



4 . 



Figure 3: The normalized longitudinal conductivity Re a xx (Cj)/eo at q e {1, 2, 3, 4} for various values 
of the flux /. The low temperature approximation (3.37) is shown as the dotted line for each q. 

We can find the permittivity e from the hydrodynamic limit T S> uj,q [38], 

_ ePuoq 2 _ STrDTCoq 2 
" Im Ctt ~ u,2 + {Dff ~ Co 2 + (vDTq>y ' (3 ' 38) 
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In this regime, the spectral function is dominated by the diffusion pole Co = —iTrDTq 2 , as 

dictated by Fick's law (A.l). The diffusion constant is D = /(/), as we calculate in 

appendix A, where we also define the function /(/). Comparing to our numerical results for 
the spectral functions as shown in figure 4 for various values of / and q <C 1, we find e = -f^j- 
We can verify the latter from the definition of the permittivity [39], 

e = 1 lim C tt , (3.39) 

which is in perfect agreement with the numerical result. 



3.2 Temperature and / dependence 

In the previous sections, we found an inter- 
esting dependence of the conductivity on 
the temperature and the flux /. These 
properties characterize the nature of the 
defect, as shown more in detail in figure 5. 
There we see that at low T or large /, there 
is a conduction threshold at to = q. We can 
interpret this as the energy required to ex- 
cite a collective excitation of the conduct- 
ing mode. In the regime u < q, the conduc- 
tivity appears exponentially suppressed as 
one might expect with a chemical potential 
q. That is, this exponential suppression in 
the low-temperature "DC limit" might be interpreted as the Boltzmann tail of some thermal 
distribution function. Examining this behaviour in more detail in the next section suggests 
the introduction of an effective temperature, which seems to play an interesting role in the 
subsequent analysis. Examining the conductivity at low T or large / also reveals "oscilla- 
tions" in the spectral curves. The frequency of this oscillations has a non-trivial dependence 
on / and seems to depend inversely on the temperature, as one might expect from the gen- 
eral scaling properties. Their amplitude is roughly independent of /, but depends on some 
positive power of the temperature and decreases with increasing In the following, we will 
also extract some quantitative approximation to this pattern. 

First, we study these two effects analytically as a perturbation around the zero tempera- 
ture limit. Next, we will show how they arise from poles in the spectral functions that can be 
interpreted in the field theory as the quasiparticle states of the resonances on the defect and 
arise on the gravity side through the quasinormal modes of the vector field. Finally, we will 
demonstrate the latter by reconstructing the location of the poles in the complex frequency 
plane from the data on the real axis and also by analytically solving a toy model that is very 
similar to our present problem. 




0.05 0.1 0.15 0.2 



Figure 4: for / e {0,4} and q e 

{0.1,0.2}. The solid lines correspond the approxi- 
mate result (3.38) with only the diffusion pole ap- 
pearing in the hydrodynamic limit. 
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Figure 5: On the left, we plot the normalized transverse conductivity Re a yy (u))/eo for q = 1 and 
/ = 2 at various temperatures, in terms of the frequency ui. On the right: Re a yy (oj)/en for q = 1 and 
T = 2 for various values of /. In both plots, the T — > limit (3.27) is shown with the narrow black 
curve. 



3.2.1 Effective temperature 

First let us study the temperature dependence of the DC limit. This can be easily done by 
finding an approximate solution in the q 3> T limit. For simplicity, we will also take u «T. 
We will concentrate on the transverse correlators, but we will see that the conductivity is 
obtained from a small perturbation around a large background, such that by (3.35), a similar 
behaviour applies to the longitudinal correlator in the limit that we will consider. Before 
proceeding, it is useful to recall here that h(u) = 1 — u 4 while H(u) is given by (3.24), so that 
asymptotically asm 0, H —> ^/^T7 I while near the horizon where u — ► 1, H —> 1. 

Let us first re-express (3.23) in terms of the Ansatz Ay = Ayoe-f such that C yy = 
, £(1 lim M ^ 

ti H' \ , H (d? 



We see that for large q (i.e., q ^ T), this equation is dominated by the second and last terms, 

h 



such that an approximate solution is £ = ±Co 5 Co = Qy jr- Implicitly, here we have chosen the 



branch corresponding to Ay decaying near the horizon. Further, as we will see below, this also 
corresponds to an infalling boundary condition at the horizon. The terms that we ignored 
are then of the order qh~ 3 ^ 2 , such that the approximation is valid in the region l-u» q~ 2 . 
The sub leading terms in ( are of the order (fifth -1 . 

Next, we study the linearized equation for a small perturbation ( — ► Co + <= : 



with the general solution 
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In the limit that we considered for Co this reduces simply to e = eoe 2 fo du <^° because of the 
exponential suppression in the last term in (3.42). As one would have physically expected, 
this perturbation grows as one approaches the horizon, and decays away near infinity. The 
subleading terms from the part of a(u) that we ignored in the integral in the exponent is 
again of order (fvPh^ 1 . 

To find eo, fix the u dependence and further constrain the subleading terms, we proceed 
by considering an approximate solution in the region h = 1 — u A <C 1, which has overlap with 
h » q~ 2 . The equation we need to solve is now 

which has a general analytic but not very illuminating solution in terms of Bessel functions, 
allowing for a combination of infalling and outgoing waves at the horizon u = 1. Choosing 
an infalling boundary condition leaves us with 



= q 2 h F 1 (2 + l f,q 2 ±) - 16,f(2 + zf) F 1 (2 + ^;g 2 A) 
8h(2 + z^) F 1 (l+zf 1 q 2 -^) 



where o^i(a;x) is the confluent hypergeometric limit function [40]. To match with the h <C 1 
regime of the asymptotic solution, we begin by expanding to first order in u> and then do an 
expansion around g 2 /i > 1, which gives us 

Hence the full solution for h 3> q~ 2 is: 

C = -q^+l+A(u,qJ) - iCoq^e- 2 " S ^(l + B(u,qJ)) , (3.46) 

where A(u,q,f) is some function that behaves away from the horizon as < 0{v?,(p) and 
B(u,q,f) behaves as < O(u ,q°). Near the horizon, i.e., for h <C q~ 2 , the solution behaves 
as C ~ —ifi — g+2ti' 2 + ' ' ' • As a consistency check in the region 1 > /i > q~ 2 , it is easy to 
verify that the (small) imaginary and (dominating) real parts do indeed satisfy (3.42) when 
taking into account the next-to-leading terms. 

From (3.46), we find that the leading term in the conductivity is 



o yy ~ 16^e ge 2 'M* . (3.47) 

Inspired by a Boltzmann factor, and by the zero-temperature conduction threshold ujq = q, 
we can interprete the exponential factor as exp[—q/T e ff] where 



(3.48) 
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We note that the integral is finite since the integrand converges as /i 1 / 2 at u — > 1. There are 
two limits in which we can evaluate this integral analytically: / = and / S> 1. In these 
limits one finds 



f = 



-eff 



-Ell) 

r(!) 



T(3/4) 2 T 
T 

7? 



1.198T 



1 + 



4r(3/4) 2 



7T' 



3/2 



(3.49) 
(3.50) 



0.847— (l + 0.372 
\l f V 



Comparing these results with the numer- 
ics, we find good convergence in a consistent 
manner of both the profile of A y {u) and the 
effective temperature with increasing q. Since 
the approximation that gave us the integrand 
in (3.48) is valid up to roughly h > q~ 2 , we 
expect that T e ff/T measured at finite q has a 
relative accuracy of roughly q~ l . A simple way 
to estimate the effective temperature from the 
conductivity is to compute dq log ~ yjj 
at large values of q. The factor of q -1 that 
we have included here ensures that the conver- 
gence to the actual value of T e f f is faster than 
logarithmic in q. In figure 6, we show the com- 
parison to the numerical estimate computed at 




Figure 6: Ratio of the "effective temperature" 
derived from the temperature and q dependence 
at q ^> 1 3> Co to the blackhole temperature as 
a function of the flux /. We show the exact ex- 
pression derived from the q — > 00 limit and the 
numerical estimate at finite q. 



q/T e ff ~ 47, that is the best numerically sta- 
ble estimate, and demonstrate how the estimates converge to the exact results. 

As illustrated in figure 6, our new effective temperature does not match the actual tem- 
perature of the system, except for / ~ 0.85. At this point, we emphasize that, as discussed 
in section 2, the degrees of freedom on the defect are in equilibrium with the thermal bath 
of adjoint fields with temperature T. Of course, T e ff is still a scale that seems to play an 
interesting role in the defect conformal field theory, as we will see in the following. Again, the 
reason that we assign this scale the appellation of "effective temperature" is that it appears 
to play the role of a temperature when the conductivity (3.47) is interpreted as a Boltzmann 
distribution. It would be interesting if one could also give a physical interpretation to the 
pre- factor 16irq in front of the exponential in (3.47). 



3.2.2 Resonances on the defect 

Next, we study the oscillatory behaviour of the spectral functions at u > q, using two different 
methods. Our results for the transverse correlator C yy obviously can also be translated to give 
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us the longitudinal correlators C xx ,Cu using (3.35). Hence we will only discuss the former 
case. 

We begin with the WKB-like expansion that gave (3.41), which was the starting point 
for the effective temperature above. Now however, we do not have the scale h ~ q 2 where we 
can match the near-horizon approximation to the asymptotic approximation. Furthermore, 

the dominant solution for e is now oscillatory since C ~ Co = ^'^' UJ \J 1 — ^|^> ra ther than 
exhibiting the exponential decay found above. The latter also reduces the validity of the 
approximation that led to (3.41) and further we have to worry about the logarithmically 
diverging integral f Co Info as u — > 1. So let us take the solution (3.42), but now with 

Co = i^fuyJl-h^, a = 2Co-(ln-^=) , (3 = £ + (in -A=) Co (3-51) 
and match this in the limit u — > 1 to the appropriate expansion of (3.44): 

< ~ "f + Sn. + ™ - -4 + + ™ ■ < 3 - 52 ' 

Now, we see that the divergent oscillations from the Q terms in (3.42) must cancel, and the 
approximation e <C Co should be valid near the horizon. Taking the limit lim M _>i e = -,^1 9 --, = 
en of (3.42), and solving for eo gives us then 

i f 1 

e\ u=0 = e = e H e- to dua ^ + / due H e~ % (3(u) 

Jo 

. w _ j\ nh ^E e -^ ( eHa _ l Co _^L) . (3 ,3) 

It turns out that the &{Yr x ) divergent terms in shcx + P in (3.53) do indeed cancel, such that 
the integral converges with the integrand oc h?K lu) / 2 as u —> 1. Unfortunately, we were not 
able to evaluate this integral analytically, even in the limits where various quantities involved 
getting large or small. We show this approximate result (3.53) compared to the full numerical 
result in figure 8. 

Because of the rapid convergence as u ^ 1, we see however that most of the contribution 
to the integral comes from regions where h ~ 1, in particular for large /. Hence as a very 
crude approximation, we can set h = 1 and hence = 0, which allows us to compute the 
integral analytically: 

eo ~ e ^vT+7e- 2 -V^A/)VwW . (3.54) 

While we do not expect this latter expression to give us the correct phase and amplitude 
information, we still anticipate that this result gives a good approximation for the frequency 
of the oscillations, 2y/l + f 2 I(f )y/l - q 2 /Co 2 - 

There is an alternative way of seeing more physically from the bulk point of view, how 
the finite temperature effects arise by casting the equation of motion for A y (3.23) in the form 
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of the Schrodinger equation, as suggested in [41]: 



(—dp+hq 2 )A y = u? A y where p 



du 



H 



(3.55) 



In terms of this new radial coordinate, the 
horizon gets mapped to p — ► oo, and we 
see that p is rapidly varying only for u < 
/ _1 / 2 and for h <C 1. This suggests that 
for large / we can approximately split the 
problem in two regions: An asymptotic one 
where h ~ 1 and p rs -' Poo (u) = Jq duy/H 
and the near-horizon region, where /i < 1 
and p ~ Ph{u) = Po — \ ln(/t) for some po- 
Going even further in our approximation, 
we assume a square potential V 
P < Poc(l 




2 

for 



q" for 



P > Poo(l) 



= VT+Pl(f) and V 
which is displayed in figure 7, 



Figure 7: The effective Schrodinger potential 
V(p) = q 2 h for the gauge field on the brane. We 
set q = 1 in the plot. 



where we see that this approximation is in- 
deed justified. At this point, we might also observe that the effective Schrodinger potential 
appearing here is very similar in structure to that found for supergravity modes [42] and for 
mesonic modes, as discussed in [43]. 

With the square potential, it is trivial to find the solution for infalling boundary conditions 
at the horizon: 



Ay ~ < 



A 2y/l - jp/tfe-tob-y/i+PHf)) 



: p < VTTpI(f) 



Keeping in mind the change of coordinates, this gives us in terms of the Ansatz that we used 
for the perturbative treatment 



C\u=0 



-VjJ 



x < 



2 (^ v / i+7M/)v / wW) 

l~q 2 /iv 2 

, cosh 2 {uy/T^pIU)yJq 2 l^-l)-qVC. 



uj < q 



The solution for u > q has the same location of the maxima as (3.54), up to a small shift 
because of the overall slope of the curve, but it is missing an exponential suppression factor 
(for increasing frequencies) in the amplitude because we approximated the smooth potential 
by a discontinuous one. For Co <C q, we also find the exponential suppression that leads to 
the effective temperature computed at u — > 0, (3.46). Hence, we can clearly see how both 
effects arise from a resonant mode on the width of the defect, and from tunnelling through 
the defect region, respectively. In appendix C, we approximate the potential by a hyperbolic 
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tangent, for which we can find an analytic solution, and find that is very closely reproduces 
the exact result with a significant deviation only at frequencies |1 — q 2 /Co 2 \ <C 1 where the 
spectral function is most sensitive to the details of the potential. 

We show a comparison of the conductivity obtained from the different approximations 
in figure 8. As expected, the perturbative approximation in (3.53) gives a very close approx- 



imation for small perturbations around the T — > result, 



<C 1, but deviates 
re analytical result 



limj^O Cyy 

significantly wherever the finite temperature effects become important. T 
(C.7) from the approximate tanh potential (C.3) however, provides a good fit for small q and 
all values of Co. For larger q > it and Co > q, there is a significant phase shift proportional to 
the separation of the resonances but their amplitudes, separation and the tailoff for Co < q fit 
very closely. This is because the phase <\> is sensitive to absolute changes in the integral of the 
potential, 8(f) oc f 5V oc q 2 , such that already small deviations in V/q 2 may have a big effect. 




Figure 8: Comparing the conductivity obtained from the numerics, from the approximation 
(3.54), from the approximation using a tanh potential (C.7) and from the sum of residues ob- 
tained from its poles (C.16). Here we focus on the oscillatory behaviour by plotting Aa yy = 
Re (& yy (oj) — £o\/l — q 2 /ui 2 ^j , i.e., we subtract off the low temperature limit (3.27). Left: / e {0,4} 
and q = 2. Right: q e {2, 4} and / = 4. 



4. Quasinormal modes and quasiparticles 

In general, the thermal correlators will have poles in the lower half of the complex frequency 
plane — e.g., see discussion in [41] or [43]. The positions of these poles characterize the 
energy and lifetime of various excitations in the system. When one of these poles is close to 
the real axis, the spectral function will show a distinct peak and the corresponding excitation 
can be interpreted as a quasiparticle. That is, the excitation satisfies Landau's criterion for 
a quasiparticle that the lifetime is much greater than the inverse energy. As illustrated in 
figures 2 and 3, which essentially plot the spectral function, the defect theory is developing 
metastable quasiparticles in the large / regime. Hence it is of interest to examine the pole 
structure of the correlators and the spectrum of quasiparticles in the defect conformal field 
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theory. This gives us not only more information on the defect field theory, but also allows us 
to speculate more on the nature of the defect. 

In principle, we could always find the poles in the thermal correlators by simply nu- 
merically computing them over the entire complex frequency plane. Of course, such a brute 
force approach would present an enormously challenging problem at a technical level. How- 
ever, since the correlators should be meromorphic, we can alternatively try to extract this 
information by fitting along the real axis, the spectral function derived from an approximate 
analytical solution of poles and positive powers - an approach similar in spirit to that fol- 
lowed in [43]. To do so, we use the complex "rest frame" frequency v = \Juj 2 — q 2 which maps 
[0,q] — > [iq,0] and [q, oo[-^ [0, oof. The motivation to do so is the fact that the resonance 
pattern is most suitably characterized by ^<Jy y — 1 as a function of u, as shown in figure 9. 
There we see that even at finite temperature this quantity varies only slowly with varying 
q. Certainly, in the low temperature limit, we expect Lorentz invariance to be restored and 
then correlators will naturally depend on the combination ub 2 — q 2 , as is implicit in (3.27) and 
(3.37). 



4.1 Finding the Ansatz 

The strategy that we will take to find the 
poles is to take a suitable Ansatz for the 
location of the n th pole in the complex fre- 
quency plane, v n = ±Dq(ti) — i^o(n), and 
for the corresponding residue, and allow 
for the parameters to vary slowly. If the 
Ansatz is good enough, and the parameters 
vary slowly enough, then we can fit the con- 
ductivity resulting from a sequence with 
constant parameters (£oj7o) to the numer- 
ical result using only the data in the region 
around the n th "resonance" . This data can 
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Figure 9: fa yy (i>) for q e {7r/4, it} and / <E {0,4}. 
This form of the spectral function is very similar for 
different values of q. 

'resonance" . This data can 
be parametrized by the amplitude of the resonance around the background and by the gap 
between the resonances. The resulting parameters (vo(n), jo(n)) then give the location of the 
pole i> n , and it's residue. 

A suitable guess for the full Ansatz is 
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v - n(z> - i7o' 



+ (^o - ilo. 



(4.1) 



where the constant terms were introduced to cancel the otherwise divergent behaviour of the 
series, since the pole terms do not decay for large n. The condition that allows us to locally 
treat the sum as an infinite series with constant {z^To} is now d n log{i>o(^)> 7o(^)} ^ 1 an d 
9„log{Po(^),7o(«)} *C ~2^i • Rewriting (4.1) in a more suggestive form, we find for PGR, 

n To 
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i.e., oj > q 

/ a \ sin IS - isinhlg% 
Im C vv = e Im z>cot vr = e Im i> , (4-2) 

"o ' '0 v ~+" '0 

such that we get the conductivity 

v sinh |p| 

ayV = £ ° & cosh l! C o S 2m ' (4 ' 3) 

which turns out to be finite at z> — > 0. These exponentially suppressed resonances are char- 
acteristically what we expect and we can, in principle, fit the parameters i>o and 7o to the 
resonance pattern. 

To be more precise however, we need to go back to the original "physical" frequency Co. 
Keeping the location of the poles and the residue fixed, the sum becomes now 

C m = e Q *tog(^£) + (4.4) 
n(D - 1%) ( n(z> - «7o) 2 n(D - ry ) 2 \ 




+ 



(q 2 + n 2 (-i) + ijo) 2 ) 1/2 \& - (<? 2 + n 2 (v> - Ho) 2 )^ 2 {q 2 + n 2 (v> - ?7o) 2 ) 1/2 , 
n(P + «7o) I n(i> + i%) 2 n(D + i%) 2 \ 



(q 2 + n 2 (D + no) 2 ) 1/2 \co + (g 2 + n 2 (D + i 7 o) 2 ) 1/2 (<? 2 + n 2 (i> + i%) 2 ) 1/2 J J 

where the term ^ log cancels the unphysical negative DC conductivity in the g> z>o, 7o 

limit that would arise otherwise. Note there is still a logarithmic divergence in the real part, 
that we are not interested in. This sequence does not sum to any known analytic expression, 
but the integral approximation can be computed straightforwardly analytically, such that in 
order to eventually study the sequence numerically, we will only sum the first few hundred 
poles and add a small "background" contribution from the rest of the poles using the integral. 
Following the same considerations, we also find an Ansatz for the longitudinal correlator, 

n 1 I {n-^)(i>o-ilo) 1 

^xx — £o 



V(g 2 + (n-i) 2 (P -i7o) 2 ) 1/2 ^-(^ 2 + (n-i) 2 (Po-i7o) 2 ) 1/2 
(n-\)(i> Q + i7o) 



{q 2 + {n-\f(^ + * 70 ) 2 ) 1/2 Co + {q 2 + (n-±) 2 (P + *7o) 2 ) 1/2 , 



(4.5) 



which converges and needs no regularization or terms with positive powers of Co. In terms of 
z>, the poles are located at (n — l/2)(±z>o — 170), with residues n _\y 2 > as we expect by (3.35) 



from the ansatz for C yy . 



In order to finally obtain the location of the poles and their residue, we split the spectral 
function at the minima into segments around each maximum and simply fit them to our 
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Ansatz giving us a set of parameters that we attribute to the local properties of the sequence 
at the most nearby pole, as described in the beginning of the section. We obtain both i> and 70 
and an overall factor (1+e^) (or (1 + e^)* on the negative branch) for the residues. The latter 
is needed because the resonance pattern is exponentially suppressed already for reasonably 
small v > it, such that the background of the fitted sequence needs to be adjusted to in precise 
agreement with the background of the data, in order to extract the relevant information which 
is contained in the resonances. As it turns out that \e&\ <C 1, we will not comment about 
its value for the rest of the paper, because it is irrelevant for both the quantitative and the 
qualitative discussions. In principle, one can introduce more parameters, such as an overall 
shift in the frequency, but this would not improve the results, since in practice, it simply 
introduces extra degeneracy in parametrizing the fit. 

4.2 Quasiparticles from the collisionless regime 

Given the Ansatz for the transverse and longitudinal correlators above, we will now discuss 
the results for determining the positions of the poles. The results are only displayed for 
Re v > but as shown in (4.1), there is a corresponding set of poles with Re v < 0. In this 
section, we focus on the collisionless or short-wavelength regime with Co S> 1 and q£,l. 

As a first test, we compare the fitted location of the poles to their exact location for the 
tanh potential in appendix C.2. We expect that this gives us a good estimate for the quality 
of the fit for the actual spectral functions. Some typical results are shown in the first plot of 
figure 10. We find that for / = 0, the fit is very poor, with the q = it/ 4 result being worse 
than the q = ir case, and there is a small deviation for / = 4 at large Co, again with a slightly 
better fit for larger q. Apart from that, i.e., for large / or large q and small Co, the fit is 
very good. This is just what we would have expected from our conditions for the validity of 
the Ansatz as smaller q imply more rapidly varying £o>7o at least for the first few poles and 
both small / and q and large Co move the poles further away from the real axis. Furthermore, 
for large Co, the amplitude of the resonance pattern becomes quickly suppressed and so it is 
subject to systematic deviations and noise. 

Now, let us look at the qualitative behaviour. We see that, as anticipated with the 
Ansatz, the poles lie roughly equally spaced on a straight line, i.e., they are resonances in a 
region of fixed width with fixed "mass" to inverse lifetime ratio. With increasing /, both the 
separation of the poles and the slope of the line of poles decreases, i.e., the poles are moving 
closer to the real axis. Of course, these changes are reflected by the appearance of distinct 
peaks in the previous plots of conductivity at large /. This behaviour is roughly independent 
of q, and there is an overall shift depending on q, that is larger for smaller values of /. One 
might expect both the decreasing energy gap and the increasing mass to width ratio since 
the length scale due to the width of the defect increases and the shape of the step in the 
potential is approximately fixed. The deviation of the poles from a straight line is stronger 
for large q and reflects the fact that the shallow potential at small q is fully probed at small 
v, whereas at large q, the resonances at small v are only sensitive to the details of the top of 
the potential, and only probe the steeper regions at higher n. This effect is obviously more 
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visible at large / because of the closer spacing of the resonances. We also see that the poles of 
the transverse correlator lie roughly half-way between the poles of the longitudinal correlator, 
as anticipated in (3.35). Comparing the various approximations to the location of the poles 
of the actual spectral function, we see the behaviour that we saw in figure 8 encoded in a 
different way. Here, we see the shift of the poles of the tanh potential that is more significant 
for larger q. 

As an aside, let us briefly return to the cuts that appeared in the the transverse and 
longitudinal conductivities, in (3.27) and (3.37), with the limit T — > 0. Given the present 
analysis, it is natural to conclude that this cut arises through an accumulation of poles near 
Co ~ q. Assuming the locations of the poles in the v plane, u n , to be roughly independent of 
q, we find that for g > 1 we get Cb n = \/v 2 + q 2 ~ q + v^/q- This then leads to an infinite 
number of poles accumulating near Co = q as T — > and resulting in a cut. 
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Figure 10: Left: Comparing the exact location of the poles of the transverse correlator to the fit, using 
the approximate tanh-potential of appendix C, at values / <G {0,4,25} and q = 1. Right: Comparing 
the poles of the transverse correlator with the ones of the longitudinal correlator at f £ {0, 4, 25} and 
q £ {1,4}. The lines are only shown to guide the eye as to which poles correspond to the same values 
of /, q- 



While the Ansatz (4.1) fixes the poles along a straight line with a fixed spacing, i.e., 
v = n(i>o— z7o) in the Re v > region, we fit the parameters locally to each peak of the spectral 
function and so the fitted poles deviate slightly from this simple Ansatz. Keeping in mind 
the limitations, let us try to extract some quantitative information on these deviations. In 
particular, at large n, the poles approach a straight line of the form v n = Su+iS^+nipQ+i^o) 
[p n = 89 + i5^f + (n — l/2)(z/fj +270) in the longitudinal case). To extract this information, we 
use different techniques in different regimes, which we outline to forewarn the reader about 
the validity of the results. For large /, where we have at least the first 5 well-fitted poles, we 
ignore the first 0-3 poles, leaving us at least 5 poles, such that we can fit the asymptotic lines 
plus a decaying exponential to Re v n and Im v n and still get information about the accuracy. 
For some cases, the exponential fit fails, and we resort to fitting a straight line and estimate 
the accuracy from the second derivative in the location of the poles. 

The assumption of an exponential deviation from a straight line may seem somewhat 
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arbitrary, but it turns out to be the right choice, as it is the only natural candidate whose 
results are independent within errors from the particular choice of the number of poles used 
for the regression. The location of the poles from the tanh potential, for example, contains 
by this criterion an (Inn) term as expected from (C.10). 

In the borderline case, where there are 4 poles, we extract the uncertainties by fitting the 
last 3 poles to the asymptotic line with the deviation estimated by the fit with 4 points. For 
3 poles only, we still get a rough estimate for the asymptotic limit (from the last 2 points) 
and for the accuracy by including the first point. For / = 0, we always find only the first two 
poles, so we can give only an order of magnitude guess for the rest of the sequence. Finally, 
we estimate the uncertainties from the errors in the fit of the sequence of poles and from the 
deviation of the estimated location of poles to their exact location in the case of the tanh 
potential. We use the latter also to add a shift to try to correct for systematic errors in the fit 
of the Ansatz (4.4, 4.5) to the numerical data. We are somewhat sloppy with the uncertainties 
in the sense that we do not distinguish between random and systematic errors. So we assume 
that the accuracy of the fits is limited by systematic uncertainties in the convergence towards 
the asymptotic straight line, which may result in a slight overestimating of uncertainties in 
the averaged data that we present below. As expected, the results from the transverse and 
longitudinal poles are identical within the errors and so we average over them. 

Let us now examine some of the results of our fitting in figure 11. In the first two plots, 
we show results for the energy gap between the quasinormal modes, £>q. In particular for 
large /, we see that the asymptotic behaviour of tt/vo matches a simple straightdine fit: 
7r/Po = ci\/y + Co with c\ ~ 1.821 and cq ~ —0.539. In the first plot, this behaviour seems 
to match well with the asymptotic behaviour of 7rT/(2T e jj) and irD(f)T. Note, however, 
that the second plot shows that upon closer examination, the deviation between tt/uq and the 
curves set by these scales in the large / regime seems to be beyond the errors expected for our 
numerical fit to t)q. Note that large / behaviour in (3.50) gives nT/(2T e ff) ~ 1.854v7 — 0.690, 
while (A. 11) yields irDT ~ 1.85Ay/J — 1. The asymptotic behaviours for these quantities have 
precisely the same slope and the difference is in the constant term (and the subleading 1/y/J 
terms), as can be seen in figure 11. This slope is only a fair match for that found in our 
straightdine fit. We expect that this is because of subleading l/y/J terms and that we would 
see better convergence at larger values of /. In any event, it seems then that Dq is closely 
related to other characteristic scales in the defect theory. Note that here since z>o appears to 
be independent of q within the errors (see figure 12), the data in figure 11 is averaged over 
q G 7r/2, 7r, 2-7T. 

We also show the overall shift £7 and the ratio 70/^0 separately in figure 11 for the cases 
q € 7r/2,7r, 2tt. In each case, these parameters show a l/y/J falloff for large /. In particular, 
this means that the width 70 is falling as 1/f and so we see the origin of the quasiparticle 
peaks in the spectral curves. In each plot, we also show irT/(2T e ff) for each case and see 
there is good agreement within the estimated errors of the numerical results. This is just 
what we expect, since the detailed shape of the step in the effective Schrodinger potential and 
hence the ratio between the two modes in the asymptotic region of the potential, is to a good 
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approximation independent of /. The slowly varying part of this ratio gives rise to the finite 
shift and the exponential suppression factor gives rise to 70/^0 > which are in the limit of large 
Co proportional to the inverse of the width of the asymptotic region of the potential. This can 
be more easily seen from the expressions in appendix C.2. From the boundary point of view, 
it comes at no surprise that the overall shift of the poles is proportional to the overall energy 
scale, and that the quasiparticle excitations become more stable with increasing /, which is 
proportional to the width of the potential step. One could make a similar plot of 59 but we 
do not show the results here. While on the whole the trends appear similar to those for 5j, 
the values are typically smaller by a factor of roughly 2 while the relative errors are larger 
by a similar factor. Hence at least for the smaller values of q, the results are consistent with 
zero shift. 
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Figure 11: Top left: The energy gap between resonances i>o presented as tt/vq. The asymptotic y/J 
behaviour is fit is several ways. Top right: The difference between the previous fits and the numerical 
results for 7r/ i> - The data in the top two plots is averaged over q £ 7r/2,7r,27r. Bottom left: The 
"mass to width" ratio 70/^0 for the sequence of poles for several values of q. Bottom right: The 
/-dependence of the overall shift of the poles Sj. The points are the average numerical data and the 
narrow lines indicate the la uncertainties, which includes both random and systematic errors. 



Now let us turn to the q dependence of the quasinormal modes. Because of the good 
agreement of the / dependence with T e ff/T, we improve the accuracy of our results by 
taking a (weighted) average over the suitably scaled values of the characteristic quantities for 
/ £ {9, 16, 25, 36}. In figure 12, we show the q dependence of the same quantities as in figure 
11. First, we see that Dq, which is supposed to depend only on the width of the potential, is 
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within the uncertainties independent of q. Any change in q however scales only the height of 
the potential step. 

Prom the results in appendix C.2, we would expect that varying q changes only the overall 
shift of the poles, but we know already that the full result has fundamentally different char- 
acteristics coming from the shape of the potential step because of the absence of a significant 
(Inn) term. In general, however, should not change significantly, since we consider here 
only / > 9, so the potential is already so wide that small details of probing the potential step 
should not change the the quasinormal modes too much. Both the shift, and the deviation 
from the linear Ansatz conspire to give us both the right "low temperature background" with 
approximately symmetric oscillations around it as in figure 9. From the fact that this be- 
haviour resembles that in the resulting conductivity from (4.2), one should assume that there 
are small shifts and deviations for small q. One also expects the shift to grow not faster than 
oc In q, provided that the ratio of the two modes in the asymptotic region depends at most 
on some power of the height of the step in the potential. 

In figure 12, we find roughly this behaviour of the shift, with small 59, 7 at small q and 
an indication of some converging or slowly growing behaviour at large q. We also find a small 
drop in 70/^0 with some converging behaviour at large q. In principle, we could try to use this 
information to try to reverse engineer the calculations in appendix C.2, i.e., to reconstruct 
the ratio of the incoming and outgoing modes. For example, the absence of a significant (Inn) 
term tells us that there is no significant v dependence, the approximately constant (in /) ratio 
shows us that the ratio of the modes has a complex phase (and also its value) but there 
is nothing really interesting to learn from this. A somewhat interesting point though is that 
the change in 70/^0 tells us that at small q, the potential "appears smoother" than at large 
Q- 

4.3 Poles in the hydrodynamic regime 

In this section, we focus on the hydrodynamic regime where uj,q 1. In this regime, the 
diffusion pole (3.38) dominates the structure of the correlators. One might wonder, why we 
have not included the diffusion pole into the sum for C yy , as in (4.1). As we will show below, 
this is because the diffusion pole disappears at a critical value of the wave-number, q c , which 
is below the values of q that we have considered to this point. Below q c there are two poles 
on the imaginary axis in the Co plane, one of them being the diffusion pole, and the other 
one at larger absolute imaginary values of aj, which decreases slowly as q grows, as shown in 
figure 13. While the diffusion pole is in perfect agreement with what we expected, the second 
pole, corresponding to rapid (i.e., on thermal scales) decay of long-range modes, is somewhat 
puzzling. In particular, it has a non-trivial / dependence at small values of q. It seems that 
for large /, the lifetime of those modes is not anymore proportional to the length scale of the 
defect, but increases less rapidly. 

At q c , there is a branch cut, and the poles move away from the imaginary axis out into the 
complex plane to turn into the first quasiparticle poles, i.e., the n = 1 poles in (4.1). Hence 
at this point, the transport changes from the collision dominated phase to the collisionless 
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phase. This is a good example of the interplay between various length scales. We can interpret 
this on the one hand as the height of the effective potential being smaller or larger than the 
inverse length scale of the defect (and hence the effective temperature) and on the other hand 
as separating between between modes smaller and larger than the size of the defect. From 
a hydrodynamic viewpoint, however this branch cut gives us approximately the mean free 
path, which is in strongly coupled systems proportional to, and of the same order as, the 
temperature scale, and we see an approximate scaling of q c with the effective temperature. 

On the right in figure 13, we compare q c with the various length scales in the problem, 
as we did before for the spacing of the quasiparticle masses in figure 11. Since we are in the 
completely opposite regime in terms of length scales of the perturbations, it is no surprise 
that there is significant disagreement between the scaling of q c and vq , but the disagreement is 
surprisingly small. In addition to the opposite limit of the size of the perturbations, the data 
in figure 11 contains only frequencies, which one can interpret as being related directly to the 
width of the defect, whereas here we consider the /-dependence of relevant values of q, which 
measure scales along the defect. Overall, it seems that in the limit of large /, g c ~ 2/(3D) or 
q c ~ 4T e ff/3. The relative factor between these two expressions is not surprising given that, 
in the previous section, we noted that D = 2/T e ff as / — > oo. Further, given our previous 
expressions for D and T e ff, we note that q c oc T/yfJ for large /, i.e., q c decreases as / grows. 
Then as the plot shows, up to an overall numerical factor, most features of the /-dependence 
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Figure 13: Left: The location of the poles on the imaginary axis below q c , rescaled by the diffusion 
constant. The dashed line is what one expects for the diffusion pole. Right: The critical values q c at 
which the purely dissipative poles disappear, multiplied with several length scales. 



of q c can be related to either of these other physical scales. 

In principle, the decreasing residue of the poles with increasing n allows us to track the 
location of the first few poles of the longitudinal correlator even further, directly by fitting 
a sequence of Lorentzians, but we will not bother about such a detailed discussion of the 
hydrodynamic regime in this paper. It is interesting to see however, how the small-g limit of 
the shift 8v + iS^f shown in figure 12 qualitatively agrees with a shift towards the bifurcation 
point. 

It is interesting to note that this pairing of the diffusion pole with a fast dissipative 
mode was also recently found in the quasinormal mode spectrum of black holes in AdS4 
[44]. However, an infinite number of pairs of poles were identified there, appearing along 
the imaginary axis. In that case, the critical wave-number at which the higher pairs meet at 
smaller and move off into the complex plane decreases for pairs higher up along the imaginary 
axis. We looked for similar higher dissipative modes in the present framework but it seems 
that the diffusion mode and its partner are the only modes appearing on the imaginary 
frequency axis. 

5. Electromagnetic duality and perturbative corrections 

At the outset of our analysis, we set M q = to maintain conformal invariance in the defect 
system. In the brane construction, this means the internal geometry is fixed and the low 
energy effective action on the effective four-dimensional brane reduces to Maxwell theory (3.7) 
with a fixed coupling (independent of the radius). 2 Hence resulting equations of motion are 
invariant under electromagnetic duality, which has interesting implications for the transport 
coefficients, as emphasized in [5]. 

2 In the D5-brane embeddings for M q / 0, the size of the internal S 2 varies and so the effective coupling of 
the Maxwell theory (3.7) depends on the radius. As explained in [5], the gauge field equations are no longer 
duality invariant and as a result the correlators discussed here are independent. 
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Given the Maxwell action (3.7), the gauge field equations can be expressed as 

VF^ = , VF^ = with % u = X - e^ pa F^ . (5.1) 

Hence we have electromagnetic duality with F^ u and satisfying the same equations of 
motion. Implicitly, we used this duality in deriving the relation between the transverse and 
longitudinal correlators (3.35), i.e., the key step was demonstrating the and Ay equations, 
(3.21) and (3.23), could be put in the same form. As in [5], this result (3.35) subsequently 
restricts the transport coefficients to satisfy 

U T (u, q) U L (Co, q) = -E 2 {u 2 - q 2 ) . (5.2) 

Since with q = 0, we have n T (u), 0) = H L (£o, 0), it follows that: 

or(a) =i n (^°) =£Q = 7T DTe or a{u) = De. (5.3) 

UJ 

That is, o~{oj) is independent of frequency and temperature. One can show that this remark- 
able result is consistent with the Einstein relation, 3 as noted already in [5]. 

However, as for any low energy action in string theory, we must expect that there are 
higher derivative interactions correcting the Maxwell action (3.7). In fact, the action (2.4) 
implicitly captures an infinite set of these stringy corrections, as would be illustrated if we 
expanded the DBI term in powers of F. This expansion would also demonstrate that these 
higher order terms are suppressed by factors of a' = £ 2 . In terms of the dual CFT, the 
contributions of these a' interactions will provide corrections to the leading supergravity 
results for a finite 't Hooft coupling. However, none of the higher order terms coming from 
the DBI action will modify the two-point correlators in the planar limit, i.e., in the large N c 
limit, because these interactions all involve higher powers of the field strength. One must keep 
in mind though that, as already alluded to in section 2.2, the DBI action does not capture all 
of the higher dimension stringy interactions. The full low-energy action includes additional 
terms involving derivatives of the gauge field strength [45, 46], as well as higher derivative 
couplings to the bulk fields, e.g., curvature terms [47, 48]. In principle, any such interaction, 
which is quadratic in F, has the potential to make finite A corrections to the correlators which 
we have studied above. 

In appendix B, we identified a particular higher derivative term which makes a quadratic 
correction (B.2) to the four-dimensional low energy action. This term makes the leading 
correction to the correlators, at least when the internal flux is nonvanishing. Including this 
term, the vector equations of motion become 

~ C(3) f 2 

V^=^L 2 VDi^, V^ = withg= 7==f - (5-4) 

27TV A y/1 + f 2 



3 See e.g. [39], section 7.4 for a suitable discussion of the Einstein relation. 
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We can recognize the higher derivative term as a string correction by recalling that L 2 /\f\ = 
I 2 . Clearly these equations are no longer invariant under the replacement: — > F^. One 
could attempt an a'-corrected electromagnetic duality by defining = (l — £L 2 n) F^ v . 
Formally treating £ as a small expansion parameter, one can rewrite (5.4) as 



vx^ = . 



(5.5) 



Hence an exchange F^ u — > X^ v does not quite leave the equations of motion invariant either, 
i.e., the sign of the £ term changes between (5.4) and (5.5). This then confirms the initial 
intuition that the a'-corrected low energy theory describing the four-dimensional dynamics 
of the vector field is no longer invariant under electromagnetic duality. Hence we can longer 
expect (5.2) and (5.3) to apply when finite A corrections are taken into account for the defect 
CFT. 

In the following, we examine in more detail the effect of this leading finite A correction. 
For simplicity, 4 we focus on the modifications to the transverse conductivity a yy . Here we 
simply present the results of our numerical calculations. The preliminary analysis determining 
the analytic form of the transverse correlator is given in appendix B. There are two distinct 
contributions to the modification of the correlator. First, since the bulk action contains an 
additional term, there are new surface terms (B.8) which must be evaluated in the holographic 
calculation. Remarkably, as described in the appendix, the net effect of this contribution is 
to shift the permittivity 

1 / 2 C(3) 



e 1- 



e 1 



2£ 



(5.6) 



v/A7T(l + / 2 )3^ 

The second modification of the correlator 
arises because the bulk equations of mo- 
tion have been corrected, as in (5.4). Hence 
the solutions for the vector are modified 
and this change of the solution alone leads 
to changes in the correlator coming from 
the leading supergravity expression (3.17). 
Now, we have some ambiguity in how we 
might define £n i n the theory with finite-A 
corrections. Recall that this quantity origi- 
nally appeared in (3.16) but above was sim- 
ply related to the conductivity (5.3) found 
in the infinite A limit. Hence a convenient 

choice, which we adopt at finite A, is: Eq = — * oo). Then our numerical results indicate 

that this second finite-A correction also shifts eo precisely as in (5.6) except for a factor of 

—3/2. The total shift is shown in figure 14 and the result seems to match precisely —1/2 

times the shift given in (5.6). 

4 Similar considerations apply to the longitudinal conductivity but the calculations are somewhat more 
involved. 
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Figure 14: The modification of eo = cr(w — ► oo) as 
a function of / from the 1/y/X corrections. 
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Given that the invariance under electromagnetic duality is lost at finite A, the frequency 
independence of the conductivity = a yy (Co,q = 0) found in (5.3) is also lost as shown in 
figure 15. Note that here we are plotting the change arising from the inclusion of the finite A 
corrections, i.e., 8a{u) = d(uj) — Eq where our subtraction includes the finite-A correction to 

eq, as described above. Note that in the figure, the factor l/£ oc is included to cancel 

the / dependence coming from the factor in front of the higher order term in (B.4). While the 
resulting conductivity shows an oscillatory behaviour, we note that the DC conductivity, i.e., 
at u = 0, is generally smaller than at high frequencies, i.e., for Co — > oo. The net difference 
is plotted in figure 16, as a function of /. As shown, the numerical results are very well 
fit with a simple analytic form proportional to / 2 /(l + f 2 ). Note that the first few points 
in this plot (including where the difference becomes positive) are not reliable, because of 
the high sensitivity to errors in 5eo, which was only computed approximately in a numerical 
calculation. 
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Figure 15: The finite-A correction to the con- Figure 16: a{u) = 0) — a(u) — > oo), the change 
ductivity 5a(uj) for various values of /. in conductivity from u> = to £j — > oo at finite 

A. 

One of the interesting features that figure 15 seems to exhibit is that the oscillations of 
a{Co) /(£,£o) for various values of / are all contained within some universal envelope, that is 
decaying with Co. In fact, this same envelope also applies for the conductivity at finite values 
of q, as illustrated in figure 17. In this figure, we are showing <5cr yy /(eo0 for / £ {1,4,9,16} 
and q G {0, 7r/4, tt/2, it, 2ir}, where 5a yy = a yy (Co,q) — d- yy (Co — > oo,q). Again, the factor l/£ 
is included in the figure to cancel the / dependence explicitly appearing in the higher order 
term in (B.4). In particular, we see here that the envelope appears to be independent of q. 
However, as shown in figure 18, for sufficiently small / (/ = 1/4 in the figure), there exists 
a critical value of q, above which the conductivity is no longer bounded by this universal 
envelope. Note that for the same values of /, the curves below the critical value of q are still 
bounded by the envelope. However, note that both £, eo oc / for large / and so the amplitude 
of oscillations in cr(w) alone is actually growing with /. 

At finite q, the leading order result (for infinite A and N c ) for a yy also exhibited similar 
damped oscillations which were confined within a certain envelope, as discussed in 3.2.2. 
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Figure 17: The extrema of -^Sa yy for q e 
{0, tt/4, tt/2, tt, 2tt} and "large""/ e {4, 9, 16}. 
We see that all the curves are approximately 
bound by some universal envelope function, 
that decays exponentially in Co. 




Figure 18: -^Sa yy for q <G {7r/4, 7r/2, tt} and 
/ e {1/4,1}.° For / = 1/4 and sufficiently 
large q, the conductivity is not bounded by 
the universal envelope (shown in thin dots). 



Comparing this previous envelope with that for Sa yy (for large / > 2), we see that the 
previous one does not depend only on Co, in contrast to the behaviour found above. The 
exponential decay of the amplitude at large Co is also slower than here than with the envelope 
for the leading infinite-A result. This would imply that the finite-A corrections become more 
and more significant at large uo, while they become less significant with increasing q. 

The "frequency" of the oscillations is approximately the same for the leading term and 
the finite-A correction. Figure 19 shows this in more detail by plotting ^5<J yy in terms of 



\Juo 2 — q 2 and also for comparison the oscillations of the infinite-A or "zero'th order" result. 
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Figure 19: Left:\/A(5CT yy in the rest-frame frequency v for q e {tt/2,-k} and / = 4, scaled to 1. The 
dotted lines show a yy — yjl — Co 2 /q 2 , also scaled to 1. Right: / = 16 



We find that there is a phase shift of between ir/2 and tt in the oscillations, implying 
that they will shift towards larger Co and decrease in amplitude. We can also see that there 
is some tendency for a smaller phase shift (i.e. less/no decrease in amplitude, less shift) as / 
and Co increase. In terms of the location of the poles, this implies a shift towards larger real 
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and imaginary frequencies and an increased spacing between the quasinormal modes, again 
more significant for large /, small q and large uj. Another point to view this is that there 
is a finite-A behaviour, that becomes important for the higher resonances. In principle one 
could quantify this more precisely by doing a perturbative treatment of the methods used in 
section 4, but we will not discuss this here. The shift eq can be absorbed into the residue. In 
terms of the potential, this implies that the potential becomes narrower, especially at small / 
and small q (or large T), which simply means that the length scale that we attributed to the 
strong coupling decreases and disappears. The uj dependence also implies that the potential 
becomes smoother at finite coupling. 




2.5 5 . 7.5 10 . 12.5 15 . f 




Figure 20: Left: log Sa yy (uj = 0)/£ (modulo a constant shift) and T j^ff) f° r 9 = 7r /2- The accuracy 
of the last two points are very sensitive to possible errors in 5e . Right: log Sa yy (uj = 0) and T q ^ 
for / = 4. log 5<j yy (uj — 0)/q 2 is shown to demonstrate the slow convergence due to polynomial factors 
of q. 

For uj — q < 0, the finite-A correction becomes quickly negative and exponentially sup- 
pressed with increasing / and q, roughly as described by the "effective temperature", such 
that the exponential suppression does not get broken but is possibly modified. We show the 
stLift &<T w^-^ as a function of / for q = and a function of q for / = I in figure 20. Recall 
that figure 16 shows the same results for q = 0. 

The form of Sa yy for / < 2 is more similar to the resonances associated with the infinite- 
A result for a yy , as the amplitude seems to decay exponentially with v and depends only 
polynomially on q, as we show in figure 18. Just as for the large /, the decay is slower 
than the one in the th order term. This demonstrates that the effects of finite A are more 
significant for small /, where the length scales are still dominated by the / = length scale, 
that we can attribute to the strong coupling. Further it shows that the length scale due to 
the "width" of the defect has a tendency to persist. 

6. Hall conductivity 

The conductivity in section 3 is diagonal reflecting the parity invariance of the defect the- 
ory. Recently AdS/CFT techniques were applied to study Hall conductivity in the three- 
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dimensional conformal field theories dual to an AdS4 background [7, 8]. The construction in 
[7] involved breaking the parity invariance by introducing a background magnetic field. Of 
course, these calculations with a dyonic black hole could be easily emulated here by introduc- 
ing additional background gauge fields on the AdS part of the probe branes — this would be 
a relatively straightforward extension of the analysis in, e.g., [49]. In [8], parity invariance is 
broken by the introduction of an auxiliary gauge field with a nonzero #-term. This construc- 
tion is closely related to the following where we produce an off-diagonal conductivity by the 
addition of a topological 0-term to the four-dimensional SYM action [50] . A related model of 
the quantum hall effect based on a probe brane construction appears in [16]. 

To introduce an xy component to the conductivity, we begin by considering the Chern- 
Simons part of the D5-brane action. In particular, the latter includes the following term: 

AI = ( 27r y N t T 5 f C l0] FAFAF (6.1) 



3! 

where C 10 ' is the RR scalar. Now the background (2.2) remains a consistent solution of the 
type IIB supergravity equations if we set this scalar to some arbitrary constant, i.e., C l0] = a. 
Of course, this choice corresponds to adding a topological 0-term to the action of the dual 
SYM theory [50]. Now if we recall the magnetic flux (2.6) on the internal S 2 , the above 
contribution (6.1) reduces to the following four-dimensional action 

AI = (2vr) 4 ^ T 5 2naq J F A F , (6.2) 

where q is the magnetic flux quantum number (2.6), indicating the number of D3-branes 
dissolved into the D5-brane. Thus upon integrating out the S 2 part of the probe brane geom- 
etry, this term (6.1) has become a topological theta-teim for the four-dimensional worldvol- 
ume gauge fields. Since it is a topological term, it does not modify the equations of motion 
(3.20-3.23) for the gauge field. However, it does produce an additional boundary term, 

AI = 2(2tt) 4 ^ T 5 aqJ d 3 a [A y (8 t A x - d x A t )] u ^ Q+ , (6.3) 

which will modify the correlators. Note that we have simplified the above expression by 
assuming that in the cases of interest (as in previous sections) the gauge fields are independent 
of y. Introducing the Fourier transform (3.13), this boundary term becomes 

AI = -i4w£ 6 s T 5 aq J d 3 k [A y (-k,u) (w A x (k,u) + k x A t (k,u))] u ^ 0+ . (6.4) 

Now following the same steps as in section 3.1, we arrive at the following off-diagonal contri- 
butions to the retarded Green's function: 

C xy = i8TT(2ir£ 2 ) 3 T 5 aquj , C ty = -i8Tr(2ir£ 2 ) 3 T 5 a q k x . (6.5) 

Note that in the T = limit, we expect this contribution to the Green's function can be 
assembled in the Lorentz invariant expression: 

AC^j, = i a e^p k p , (6.6) 
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where a is the dimensionless constant: 



a = 8ir(27r£ 2 s fT 5 aq = 



2aq 8aq 



(6.7) 



*9. 5ym ' 



The corresponding analysis in the D7 framework gives 



a 7 = 16ir 2 (27r£ 2 s ) A T 7 aq 7 = 



2a q 7 8 a q 7 



(6.8) 



While in principle, this form (6.6) need not be preserved at finite temperature, our results (6.5) 
calculated at finite T indicates that in fact the form is preserved. Of course, this independence 
of the temperature is undoubtedly related to the topological nature of the #-term which is 
responsible for this off-diagonal contribution. It is amusing to note that since q (and q 7 ) is an 

2 

integer, (6.7) and (6.8) take just the form of the integer quantum Hall effect, i.e., a xy = 
with n £ Z. By this analogy, we would associate e 2 = 4a/<7 s . 

7. Discussion 

In this paper, we used holographic techniques to investigate the transport properties of certain 
defect CFT's. In particular, we studied matter on a (2 + l)-dimensional defect emersed in a 
heat bath of (3 + l)-dimensional N = 4 super- Yang-Mills plasma. Our analysis covers two 
distinct defect CFTs. The first was realized by embedding N f probe D5-branes in AdSs x 
S 5 , as described in section 2.1 and in this case, the system (at T = 0) preserves eight 
supersymmetries. The second system involves embedding N f probe D7-branes in the AdSs x S 5 
background and the resulting defect CFT preserves no supersymmetries. In both cases, the 
theory could be deformed by introducing an additional internal flux on the probe branes. In 
the dual CFT, the defect then separated regions where the rank of the SYM gauge group 
was different. As described in section 2.2, this flux was crucial to remove an instability which 
would otherwise appear with the D7-brane construction. Perhaps surprisingly, the transport 
properties of both defect CFT's were essentially identical. 

Overall, our analysis revealed the expected diffusion-dominated hydrodynamic limit at 
small wave-numbers and we found a smooth crossover to a collisionless regime at the large 
wave-numbers. In the latter regime, the defect theory exhibits a conduction threshold, given 
by the wave-number q of the current and the system is approximately described only in terms 
of the "rest-frame frequency" v = \J u) 2 — q 2 . 

In many respects, our results coincided with those in [5], where holographic techniques 
were used to study a purely (2 + l)-dimensional system with sixteen supersymmetries. Hence 
maximal supersymmetry (or supersymmetry, in general) does not seem to be a key feature 
for producing the interesting behaviour of these holographic models. Instead many proper- 
ties seem to emerge from the infinite iV c and infinite A limits, that are implicit in making 
a supergravity analysis of the AdS dual. In section 5, we elucidated one such effect that 
arises purely from the large- iV c and large-A limits, namely the frequency independence of the 
conductivity, a(oj). In [5], this effect was described as a consequence of the electromagnetic 
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duality of the gauge theory giving the dual description of the CFT currents. We were able 
to explicitly show that this duality is lost when stringy corrections are included in the world- 
volume gauge field action and explicitly calculated the frequency dependence in a (to) arising 
from the corresponding finite- A corrections to the conductivity. As described in appendix B, 
one can well imagine that there will be other interactions which, although they appear to be 
of higher order in the a' expansion, provide further corrections to the conductivity f n /y/\, 
where n > 2. Hence, our results in section 5 are only the leading corrections when / is small 
but finite. There will also be curvature interactions to the worldvolume action of the probe 
branes [47, 48]. These will also produce fmite-A corrections but in contrast to the previous 
discussion, the latter will not be enhanced by factors of / and first appear only at order 1/A. 
While a completely consistent set of curvature terms is not known at this time, the effect of 
certain representative terms was considered in [24]. 

While certain aspects of charge transport were similar for the present defect CFT's and 
the maximally supersymmetric CFT studied in [5], we also found some profound differences. 
The most prominent is the dependence of our results on the internal flux /, certainly a dif- 
ference since no such parameter exists in the maximally supersymmetric case. For example 
then, with a large /, strong quasiparticle peaks appeared in the spectral functions and con- 
ductivities. Similarly, certain phenomena in the defect CFT seemed to be controlled by a new 
dynamically generated scale in this regime, i.e., a scale distinct from the temperature T. We 
denoted this scale as the effective temperature T e ff in section 3.2.1. For small /, T e jf ~ T 
to within factors of order one. However, for large /, one finds that irT / (2T e f f) ~ k \f] where 
k = 4r(5/4) 2 /y / 7r, as shown in (3.50). While this seems a curious way to present the ratio 
of T and T e ff, it was found in section 4.2 that precisely this combination played a role in 
fixing the spacing of the quasiparticle poles. Further, as also noted there, precisely the same 
behaviour was found at large / for the diffusion constant: tt DT ~ k yJJ. 

A more intuitive picture as to the origin of this dynamical scale comes from considering 
the probe brane geometry, as in section 2. Recall that when / is nonvanishing, the defect 
separates a region where the SYM has gauge group U(N C + q) from one where the gauge 
group is U(N C ). However, on the side where the rank of the gauge group is enhanced, the 
defect also excites a noncommutative configuration of adjoint scalars in a U(q) subgroup of 
the full U(N C + q). At T = 0, this configuration extends to infinity with Tr(<£ 2 ) = ^ 4y. 
In particular, this configuration preserves the conformal symmetry, i.e., does not introduce a 
new scale. The scalar profile is reflected in the radial profile of the probe D-branes which also 
extends out to z — > oo when T = and / is nonvanishing. However, when the temperature 
is nonvanishing, the probe branes fall into the horizon at a finite value of z. For large /, one 
finds 

ZmaxT = k^ff, (7.1) 

where k is precisely the same constant found above. The natural interpretation of this profile 
is that at finite temperature, thermal fluctuations decohere the scalar fields at some distance 
away from the defect. That is, at finite T, the defect can only excite a coherent configuration 
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of scalars out to z max . However, at small /, z max vanishes whilst the dynamical length scale 
doesn't vanish. We can interpret that by some extra contribution to the width that arises 
from excitations in the bulk fields that are induced at strong coupling from the presence of 
the charged matter degrees of freedom on the defect. 

One interpretation then is that the defect effectively acquires a finite width when T is 
nonvanishing. This intuitive picture may seem more reasonable, if we recall that the system 
is at (extremely) strong coupling and so any probe exciting of the defect fields will actually 
excite a complicated mixture of defect fields and "bulk" SYM fields. This picture of finite- 
width for the defect seems to match well with the results for the quasiparticle spectrum on 
the defect. In particular, we found that both the conduction threshold and the resonance 
peaks are well-described by a quasiparticle "tower" with equally spaced "rest frame" energy 
and constant "mass to inverse lifetime ratio". The length scale that is suggested from this 
spectrum is very similar to itDT and irT / (2T e f f) (and at large / also similar to z max ) plus a 
small constant. 

In terms of the effective Schrodinger analysis, e.g., see appendix C, the quasinormal 
spectrum arises in the gravity side from interference on a potential barrier in the asymptotic 
region. From the profile of the brane, this translates into an interference or resonance in 
the region around the defect. These two dual pictures of the origin of the spectrum seem 
orthogonal. The intuitive CFT picture involves a width and implicitly, a potential, in the 
space transverse to the defect, while the effective Schrodinger analysis constructs an effective 
potential in the radial or "energy scale" direction. It would certainly be interesting to have 
a clearer connection between these two descriptions. 

The hydrodynamic and collisionless regimes are cleanly separated at a critical wave- 
number where the diffusion pole disappears, as observed in section 4.3. There, we found 
that the diffusion mode is partnered with another dissipative mode on the imaginary axis. 
These two poles coalesce at the critical wave-number and move out into the complex plane for 
larger q. Hence precisely at the critical wave- number q c , the corresponding thermal correlator 
will exhibit a curious double pole on the imaginary axis. Interestingly, q c has a similar 
qualitative / dependence as the other (inverse) length scales and is also quantitatively close, 
as q c /(irT) ~ 0.67 '/ \k^ff). This supports again the concept that the properties of the defect 
are controlled by one fundamental length scale, that can be related to the effective width of 
the defect. 

In section 6, we outlined a topological Hall conductivity, which is induced when the defect 
is coupled to the SYM gauge theory with a topological 0-term. Of course, it would also be 
interesting to study the Hall conductivity induced by a background magnetic field on the 
defect. Other interesting directions would be to study the effect of a finite mass for the defect 
fields or of a finite chemical potential. With these additional parameters, there are many 
non-trivial physical properties of the defect related to condensed matter physics that one 
can study through the transport properties. We will address these questions in an upcoming 
paper [51]. 
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A. Diffusion constant on the defect 

The worldvolume gauge field corresponds to a conserved current on the defect in the dual 
CFT. In the hydrodynamic regime, one then expects to see the diffusion of the conserved 
charge according to Fick's law: 



This expectation can be confirmed in a holographic context [52, 53, 54] and, in fact, the 
computation of the diffusion constant D can be performed in a number of different ways. In 
the following, we use the membrane paradigm approach. 

The computation of the diffusion constant via the membrane paradigm was discussed in 
[52] where explicit formulae for various transport coefficients in terms of metric components 
for a wide class of backgrounds were derived. There, the authors considered perturbations 
of a black brane background and a formula for the diffusion constant (eq. (2.27) in [52]) 
resulted from a derivation of Fick's law. An analogous computation can be performed for the 
D5-branes' vector field for black hole embeddings considered here, with the result 



In the first expression above, the metric g is the induced metric on the D5-branes (2.7) and ^mt 
is the determinant of the metric on the internal two-sphere (with unit radius). The integral 
can be evaluated analytically yielding a hypergeometric function: 



Figure 21 shows a plot of /(/). 

This same integral in eq. (A.2) reappears at various points in our analysis and so it is 
useful to gain some better intuition for this expression. First, let us rewrite the integral as 
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Now, we find that we can expand the integrand around s = as 



and around s = oo 
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At s = 1, the convergence of both sequences goes as ^ and a more precise approximation is 
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(A.7) 



Combining (A. 4) with these expansions, we can find the integral in various approximations 
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where c\ = 



ds 



VT+~s 



27^/^ (f) and Coo =/ c 



v / I+? r 



2ci ~ 1.854. Further, we note 



that the expansion about s = is just the expansion of the hypergeometric function as a 
hypergeometric series. 

Hence, at small and very large /, we 
have for the diffusion constant 
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We note that the / = result is different 

from but close to the value found for M2- 

brane hydrodynamics [55]: DT = 3/4tt. 

Also note that in the limit of large /, the 

diffusion constant grows as Jl . „. ni . , Tl r , , ,, , r 

v J Figure 21: I he integral plotted as a function 

of \ff: The solid line is the exact result as given 
B. Corrections to D5-brane action in eq. (A.3). The upper dashed (lower dotted) line 

corresponds to the first term (first two terms) in the 
The worldvolume action (2.4) is a low-energy i arge / expansion in eq. (A.9). 
effective action which captures the interac- 
tions of the massless open string modes supported on the D5-branes. The "stringy" nature 
of the underlying theory will in principle produce an infinite series of higher dimension terms 
that are suppressed at low energies by the inverse string tension, i.e., a' = However, in 
practice, one typically only includes a specific set of terms to a certain order in the a' expan- 
sion. In fact, the action (2.4) implicitly captures an infinite set of these stringy corrections, 
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as can be seen by the explicit factor of I 2 accompanying the gauge field strength in the DBI 
action. Further, we might add that this entire series of terms plays a role in our analysis, as 
the full square-root form of the action is used in (2.8) to normalize the DBI contributions. 
However, as already alluded to in section 2.2, the DBI action does not capture all of the higher 
dimension stringy interactions. Beyond the nonabelian commutator terms referred to there, 
the full low-energy action includes additional terms involving derivatives of the gauge field 
strength [45, 46], as well as higher derivative couplings to the bulk fields [47, 48]. In terms 
of the dual CFT, the contributions of these a' interactions will provide finite A corrections to 
the leading supergravity results. 

In the following and in section 5, we focus our attention on a particular new term involving 
derivatives of the field strength, which modifies the vector correlators on the D5-branes. Our 
results below give the leading 1/A corrections when / is finite. Our calculations consider 
explicitly those leading contributions for the transverse correlator. As we discuss in section 
5, the higher derivative interaction also upsets the electromagnetic duality on the AdS4 part 
of the worldvolume. 

An effective action to describe open string gauge fields at higher order in a' has been 
extensively studied in the literature, e.g., see [45, 46]. However, since we consider only the 
linearized equations of motion (or the quadratic action) to compute the correlators in the 
AdSi directions, the leading term at (a' 3 ) is [45] 



In fact, this term only becomes relevant for the calculation of the vector correlators because 
of the flux background on the compact S 2 . After integrating over the internal two-sphere, 
this interaction reduces to 



Note that the background flux (2.6) remains unmodified by this new interaction. When con- 
sidering linearized fluctuations of the vector field, we can work with a perturbative expansion 
in £ 2 /L 2 = 1/y/X. We begin with the Ansatz 



The equations determining the leading order field are still the same Maxwell equations 
(3.10), while for i^ 1 ), we have 
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Writing this out explicitly for fi = y gives 



_ c(3)/ 2 , att /» ft , /r 4(0) , fl 



Note that we adopt the convention above that the derivatives d u act on all factors to their 
right. At the horizon, u — > 1, we again wish to impose infalling boundary conditions. If we 
substitute the expansion (3.28) for A y as u — > 1, the right hand side of (B.5) diverges as /i la; / 4-2 
and hence we expect to find singular behaviour in A^ there. Hence we begin by isolating this 
singular behaviour in a particular solution to (B.5): Ay ] ' NHG ~ h^/ 4 (±+b+c\ogh +dh log/i) 
with appropriate constants a,b,c, d. This particular solution is holds to order h° near u — > 1 
and is well-behaved in the rest of the geometry. Implicitly, it also satisfies the desired infalling 
boundary conditions. Next we add to A y a contribution which is regular at the horizon 

and takes the form on an infalling homogenous solution near u — ► 1 as described by (3.28) 
and (3.29). That is, we use the Ansatz 

4," = 4 > -" HG + ^ (i+(fir£ + irk) (l - »>) #( ") ■ (R6) 

where d u ^{u)\ u ->i = 0. This Ansatz is constructed so that J^{u) is well behaved everywhere 
and we proceed by calculating this profile numerically. In practise, we apply the boundary 
condition at some small, but finite, (1 — u), so need we go one order higher in h and increase 
the accuracy for solving A°) . 

The correlator is then found by substituting our original Ansatz (B.3) into (3.17), 



e d u (Ap+ei/L*AV) ^ so 
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(B.7) 

From this expression, we can see that the normalization of any homogenous solution in A^ 
does not effect the correlator, as we expect for a gauge- invariant quantity. The expressions 
for Cu and C xx following from (3.18) are similar to that above. 

The second term on the right-hand side of (B.7) yields the correction to the correlator 
due to the fact that the solutions for the fluctuations A^u) are modified. In addition, (B.2) 
also contributes to the overall value of the bulk action and hence provides an additional 
modification of the correlators. Proceeding as in (3.12), we now get two boundary terms at 
order since D^F a pD^F al3 = D^{F a pD»F a P) - D a {ApD^F a f 3 ) + A^e.o.m.f , where 
the last term combines with contributions from the leading Maxwell term in the action to 
vanish by the equations of motion. These expressions lead to a number of new contributions 
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to the flux (3.16), which in principle even contribute to off-diagonal correlators. We might add 
that there is a further ambiguity in these expressions because the effective action (B.2) was 
constructed from examining string scattering amplitudes [45] and so it is only determined up 
to total derivatives or boundary terms. Explicitly comparing [45] and [46], one finds that in 
fact their results differ by such a boundary term. However, this ambiguity does not contribute 
in our background (3.11) and in fact, of the myriad of potential boundary contributions, only 
a single term survives 



AC il fC(3) e 



V 2 d u A y 



(B.8) 



Similar expressions survive for Ctt and C xx while the potential contributions to the off- 
diagonal correlators vanish. Combining the results in (B.7) and (B.8), we can write the 
total change to the correlator as 



y ~ l 2 ^TTP 



d u 



A 



(i)' 



(0) 



L a / 2 C(3) V 2 d u A^ 



4°y 2VIT72 4 



(B.9) 



Again similar expressions arise for 5Cu and SC XX . 

The two previous formulae, (B.8) and (B.9), still require a precise definition for V 2 d u A y . 
This expression should understood as the covariant tensor expression V 2 F uy which when 
evaluated asymptotically yields a remarkably simple result: 

(B.10) 

that they will not contribute to the 



V 2 d u A y — 1 + p d^Ay H 



where the implicit terms decay rapidly enough as u 
correlator. Hence (B.8) can be greatly simplified to 



AG 



/ 2 C(3) 



£0 



d u Ay 



mi 



£ 2 vr(l + / 2 ) 3/2 V / l+7 f A v 



(B.11) 



The last factor has exactly the same form as the leading correlator (3.17) and so this contri- 
bution can be interpreted in terms of a rescaling of the pre-factor £q: 



£0 



eo 1 + 



1 



/ 2 C(3) 



V^vra + J 2 ) 3 / 2 



(B.12) 



where we have replaced £ 2 /L 2 = l/\/X Remarkably our numerical calculations show that 
the first contribution to 8C yy in (B.9) also produces a shift of £0 with precisely the same / 
dependence - see figure 14. 

At this point, several comments are in order. We already pointed out that the square- 
root form of the DBI action already incorporates an infinite set of stringy a' corrections. 
While this form was incorporated in our leading order calculations, e.g., (2.8), it did not 
appear to introduce any 1/A corrections. Of course, these factors are hidden in the definition 
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(2.9) of / and for finite values of /, we are actually introducing a magnetic flux quantum 
number q ~ 0(y/X). In this context, it is not quite correct to say the interaction (B.l) is 
the leading term to modify the correlators. One can well imagine that there will be other 
interactions which, although they appear to be of higher order in the a 1 expansion, will modify 
the correlator with contributions of order f n /\f\ where n > 2. Of course, these contributions 
will be suppressed in a regime where / < 1, Considering the possible tensor structure of 
the relevant higher order interactions, it seems that this class of contributions will always be 
arising from an equation of motion of the form (B.4). 



C. Hyperbolic tangent potential 

To get more insight into the appearance of the finite temperature effects in the frequency 
dependence of the conductivity, we can study a qualitatively similar problem that has an 
analytic solution. 

C.l Finding the spectral curves 

Let us modify the effective Schrodinger equation in (3.55) to write it in terms of the complex 
frequency v of section 4, i.e., v = \/ uj 2 — q 2 , 

{-d 2 + V)A y = D 2 A y (C.l) 

with V = —q 2 u 4 and p = / du 

Jo 



h 



Recall in terms of the new radial coordinate p, the asymptotic boundary is mapped to p = 
and the horizon, to p — > oo. In fact, it is straightforward to evaluate the integral above to 
find p(u) in terms of incomplete elliptic integrals of the third kind or alternatively, in terms 
of the Appell hyper geometric function [56] p = \/l + f 2 u F\ Q; |, 1; |; —f 2 u 4 , u 4 ) . Note that 
the only difference between (3.55) and (C.l) is that we have subtracted q 2 A y from both sides 
in the equation above. Hence in the present form, the effective potential V vanishes at p = 
and approaches — q 2 as p —> oo. 

The equation of motion in the form of (C.l) can be examined in three distinct regions 
which for / » 1, where the potential looks roughly like 



V ~ 



V ~ -f(2f/ 2 -p)- 4 

V „ _g2 + g ~2 e -4(p-2/V 2+1) 



1/2 



f 1/2 <p< 2/ 1/2 - 1 

2/ 1/2 -1<P- (C2) 



Of course, the full potential is smooth across these three regions. 

To model the "smooth" step with a finite slope and a fast tailoff in the asymptotic regions, 
we will use 

v _ -2 1 + tanh n(p - p ) 

Vapp — —Q ^ ; V^-^J 
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which allows us to find exact solution for (C.l). We fix the "step" position po either as 
location where the the potential has half its minimum, i.e., at u\ al j = 1/2 or the point where 
the slope of the potential is minimum, i.e., at 

_ 5/ 2 - 7 + ^(7 - 5/ 2 ) 2 + 108/ 2 



18/ 2 



(C.4) 



To fix the parameter n that we will use to characterize the slope of the step of the poten- 
tial, we simply take the slope of the full potential at either of the corresponding points: 



"2 dpVL 

" Hhal f ,min 



u d h(u)H(u) 



-1/21 



> u half, 



Working with Uhaif or u m i n will give 



V app which is a good approximation of the full potential V. As shown in figure 22, the poten- 
tials constructed with these two choices typically can not be distinguished. For comparison 
purposes, the figure also indicates two other interesting scales: p = irT/2T e ff and irDT. 
There we see that both of these scales are close to the width of the effective potential but 
p = tt D T is particularly close to the center of the step. 

Note that V(p = 0) = and d p V(p = 
0) = while neither of these properties holds 
for the approximation (C.3). We will take 
account of this fact simply by expanding 
around p — po — > — oo and discarding higher 
order terms when computing the values at 
p = 0. One can expect this potential to 
be a good approximation for v S> 1, where 
the result is less sensitive to the details of 
the potential away from it's maximum slope 
and for large /, where the results should be 
dominated by the large flat part of the po- 
tential. 

Using infalling boundary conditions at 
p — > oo, the solution is 




Figure 22: The exact potential V (solid grey) and 
the tanh approximation V app (black dotted), for / € 
{0,4,25} with q = 1. We also indicate a number 
of interesting values of p with narrow vertical lines: 
p{uhalf), solid black; p{u min ), solid grey; p = 
dashed black; and p = tt DT, grey dashed. 



2T„ 



A y0 2Fl ^l + i£(l-^V^), § 1+e2 n( P -P0) 

(1 - tanh n(p - Po ))~f V^ 2 ^ 2 (1 + tanh n ( p _ po ))f 
and the asymptotic expansion gives us 

'iTW^-i) r(i(i-Vi + 9 2 /^ 2 )) 2 r(-f) 



(C.5) 



Ay — Ay 



— A,,, 



Ju(p-po) 



\ 



(yrrw^+i) r(-g (y/r+wi^+i)) r(f) 



+ e 



2iu(p-p () ) 



(C.6) 



for which we redefined A v q. In the opposite limit, as p — > 00, the solution converges expo- 
nentially to A y oc e~ wp v /l +« 2 / i/2 , which is the desired infalling wave solution at the horizon. 
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The transverse conductivity can now be written in a compact suggestive form in terms 
of the implicitly defined N: 

y N _ g2wpo 

~a yy = £oRe ^ K + e2 .. po , (C.7) 

which is of the same form as the result (3.2.2), but with different coefficients. We again see how 
the oscillatory pattern arises from interference in the asymptotic region due to the potential 
"step" , and how the effective temperature arises from tunnelling through the potential out of 
the "heat bath" in the near horizon region. For \v\ 3> n, we find the limit 

N |i> ^° ^e s s n ( Re + log v/xi)) , (C.8) 

sgn(Re v) irq z 

This result shows that the subleading term in the conductivity (beyond the low temperature 
limit (3.27)) has an exponential frequency dependence for Co > q, as we found in section 3.2. f. 
We can also see how the results depend on n, i.e., on the slope of the step in the potential. 
We compare the present results (C.7) for the conductivity to our numerical results and to the 
perturbative approximation in figure 8. There we can see that the various approaches are in 
very close agreement for Co > q. 

C.2 Exact pole structure 

We can look at the solutions (C.6) and identify the quasinormal modes by imposing the 
asymptotic boundary condition: A y (p = 0) = 0. The quasinormal frequencies are then 
simply determined by the equation 

= (v/TTIW-i) r(j(i - ^Pl)' r(-f ) ^ 

(yrrsw-i) r(-i(i + v / r+?75> 5 )) 2 r(f) 



which simplifies in the limit of S> n, q to 



; (2ip -sgn(Re i>)7r/n)i> _ 4nP _ 

sgn(Re v) irq 2 



(CIO) 



These quasinormal frequencies also give the location of the poles in the spectral function 
[36, 57], with the exception of the asymptotically constant solution, v = 0. In the terminology 
of section 4, the present approximation yields: 

u o = tit - ; — m? ' ~ = n — • (C- 11 ) 
4pg + vr 2 /n 2 v 2p n 

Further as mentioned in section 4, the subleading contributions are logarithmic in n, giving 
us 

x - 7o-^0/, . , 4n(i7o-P )N mm 

bv n = — (lnn + ln ^ ), (C.12) 

2tt nq z 

in contrast to the constant shift as the largest subleading term for the full correlator. 
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The residues in the Green's function, before taking the imaginary part, are given by 



ag> ee Res 



(C.13) 



u ^o+ 2 ^o - dy log 7 1 y =i , n ' 
One finds that da log 7 can be expressed in terms of digamma functions, and asymptotes to 



£0 



Res 



d u A % 



y/T+f2i>=Un A y 



-2v„ 



~ e 



2po + ivr/n ' 



(C.14) 



where we used implicitly n G Z to label all the poles. 

We can use the knowledge of the poles to verify that the spectral function is indeed 
approximated extremely well by the regularized sum of poles plus the term linear in uj as 
given in (4.4), 



a 



VI) 



. £^ nSgn(Re "' 



1 



1/2 



,V2 



+ 



sgn(Re u n 



1/2 



+ Eq lim — log 



(C.15) 



or from (4.5) 



C X x — £o ^2 \ 



n ) Vn sgn(Re P„ 



(q 2 + P 2 ) 1/2 £ _ S gn(Re v n ) (q 2 + v 2 



1/2 



(C.16) 



which we show in figure 8. We also use the exact location of the poles to check how well the 
algorithm from section 4 finds the location and residue of the poles, which we show in figure 
10 and discuss in section 4. 
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